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ABSTRACT 

Most astronomers now accept that stars more massive than about 9 M Q explode as supernovae and 
leave stellar remnants, either neutron stars or black holes, with neutron stars being more prevalent. 
Recent modeling of the explosions suggests a significant diversity in the key natal properties — rotation 
rate, velocity, and magnetic field strength — of the resulting neutron stars that account for the association 
of active radio pulsars, pulsar wind nebulae, and magnctars with supernova remnants (SNRs). The 
discovery of a central X-ray source in Cas A, the youngest known Galactic SNR, dramatized the expected 
diversity. However, less than half of the SNRs within 5 kpc have identified central sources, and only 
three are identified as the remnants of Type la SNe. Here, we report a systematic effort to search 
for compact central sources in the remaining 23 SNRs of this distance limited sample. Our search 
was inspired, on empirical considerations, by the enigmatic faint X-ray source in Cas A; motivated, on 
theoretical grounds, by the expectation that young neutron stars emit cooling X-ray emission; and made 
possible by the superb angular resolution offered by the Chandra X-ray mission and the sensitivity of 
the XMM-Newton mission. 

In this first paper we report Chandra observations of four SNRs (G093.3+6.9, G315.4-2.3, G084.2+0.8, 
and G127.1+0.5). We have undertaken a systematic optical/IR identification program of the X-ray 
sources detected in the field of each SNR. Foreground (flare stars, active stars) and background (active 
galactic nuclei) sources have identifiable IR/optical counterparts. In contrast, the counterparts of neutron 
stars (or black holes) are expected to be very faint. We are able to account for all the well-detected 
X-ray sources and thus able to state with some confidence that there are no associated central sources 
down to a level of one tenth of that of the Cas A central source, Lx < 10 31 ergs s _1 . We compare 
our limits with cooling curves for neutron stars and find that any putative neutron stars in these SNRs 
must be cooling faster than expected for traditional 1.35 Mq neutron stars and that any putative pulsar 
must have low spin-down luminosities E < 10 ergs s _1 . However, our limits are unable to constrain 
the presence or absence of more unusual options, such as relatively more massive neutron stars with 
M > 1.45 M , neutron stars with exotic interiors, or quiescent black holes. In subsequent papers, we 
will report on the X-ray and optical/IR observations of the remaining members of the 5-kpc sample. 

Subject headings: pulsars: general — stars: neutron — supernova remnants — X-rays: stars 



1. INTRODUCTION 

Understanding the deaths of massive stars is one of the 
frontiers of modern astrophysics. Considerable observa- 
tional evidence substantiates the idea that stars below 
M w ~ 8Mq end their lives as white dwarfs (Weidemann 
1987), while the detection of neutrinos from SN 1987A 
dramatically illustrated that more massive stars undergo 
core collapse (Hirata et al. 1987). The outcome of core col- 
lapse can either be a neutron star or a black hole (Woosley, 
Heger, & Weaver 2002). However, there are great uncer- 
tainties in the mapping between initial mass of the star and 
the end product, and even more uncertainties in the natal 
properties of the stellar remnant. It is these uncertainties 
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that give observers opportunities to make new discoveries 
and theorists to predict or "postdict" these discoveries. 

The first issue — the state of the star prior to collapse — 
is very dependent on the mass loss history of stars (Heger 
et al. 2003), a phenomenon that is poorly understood and 
can easily be modified by the presence of a binary com- 
panion. For solar metallicity (the situation relevant to this 
paper), stars with between Mi ~ 9M Q and M u ~ 25 Mq 
are expected to form a neutron star (Heger et al. 2003), 
while stars above M u are expected to form a black hole 
either by fall-back of material which transmutes the neu- 
tron star to a black hole or by direct collapse. As an aside 
we note that progenitors with masses between M w and 
Mi may form an O-Ne-Mg white dwarf that may collapse 
to neutron stars or simply explode (Miyaji et al. 1980; 
Nomoto 1984, 1987). 

The second issue — the natal properties of the stel- 
lar remnant — apparently involves delicate physics but 
has strong observational ramifications. The gravitational 
binding energy of a neutron star is 10 53 erg, of which only 
1% appears to be coupled to the ejecta (which ultimately 
powers the SNR). Even more minuscule fractions go into 
rotational energy, kinetic energy (bulk motion) and mag- 
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netic fields. It is now generally agreed that three dimen- 
sional effects in the explosion determine the natal proper- 
ties (Burrows 2000; Kifonidis et al. 2003). 

So far the discussion has assumed that the only gross 
parameter of interest is the mass of the star. However, it 
is likely that rotation of the progenitor can also profoundly 
affect the outcome. 

These two issues are now propelling two different ar- 
eas of inquiry. The relationship between the progenitor 
properties (mass, rotation) and the gross outcome of core 
collapse (neutron star or black hole) is observationally be- 
ing determined by systematic studies of supernovae and 
GRBs and their interconnection. The second area is in 
understanding the natal properties of neutron stars, which 
is the main focus of the paper. 

The discovery of pulsars in the Vela SNR (Large, 
Vaughan, & Mills 1968) and the Crab Nebula (Staelin 
& Reifenstein 1968) made concrete the suggestion that 
core collapse results in neutron stars (Baade & Zwicky 
1934), some of which manifest themselves as radio pul- 
sars. Young pulsars, in addition to pulsing in the radio, 
can and usually do power synchrotron nebulae (Weiler 
& Panagia 1978) that are indirect markers of pulsars. 
These synchrotron nebulae are commonly called pulsar 
wind nebulae, or PWNe. Over the following two decades, 
the notion that neutron stars resemble the Crab pulsar 
guided the search for central objects as well as intensive 
radio mapping of SNRs. As a result of these efforts, the 
term "composite" SNR (PWN + shell) was added to the 
SNR lexicon (see Milne et al. 1979). 

However, over the last five years there have been three 
developments that have severely revised our picture of 
young neutron stars. First, astronomers have come to 
accept of tremendous diversity in the natal properties of 
young neutron stars. Anomalous X-ray pulsars (AXPs), 
soft 7-ray repeaters (SGRs), nearby thermal and radio 
quiet neutron stars, long period radio pulsars with high 
inferred magnetic fields (HBPSR) are now routinely found 
in astronomical literature. These new classes of neutron 
stars have primarily come from high energy (X-ray and 
7-ray) observations. Second, there is appreciation that 
the radio luminosities of PWNe is poorly dependent on 
the spin-down luminosity of their central pulsars. For 
example, energetic pulsars may have faint PWNe (e.g., 
PSR J1119-6127; Gonzalez & Safi-Harb 2003), and very 
sensitive imaging of the regions around identified cen- 
tral sources has frequently yielded only upper limits to 
the surface brightness of putative PWNe (e.g., Gaensler, 
Bock, & Stappers 2000). Third, seven nearby cooling NSs 
have been identified (Haberl 2003) through ROSAT ob- 
servations. Many of these neutron stars do not appear 
to be evolved versions of standard radio pulsars; e.g., 
RX J0720.4— 3125 has a period longer than almost any 
known radio pulsar but has a typical i?-field strength (Ka- 
plan et al. 2002a; Zane et al. 2002). 

While this diversity is clearly demonstrated observation- 
ally, theory and simulation cannot yet constrain the fun- 
damental birth properties of neutron stars (e.g., Burrows 
et al. 2003). Models still have difficulties achieving explo- 
sions, much less following the activity in the post-collapse 
object in any detail. 

Three years ago we began a program of observation- 



ally investigating the stellar remnants in a volume-limited 
census of Galactic supernova remnants. The approach we 
took was inspired by the first light picture of Chandra, the 
discovery of a central X-ray source in the well-studied and 
youngest known supernova remnant in our Galaxy, Cas- 
siopeia A (Tananbaum 1999). The nature of the object 
continues to be debated (Pavlov et al. 2000; Chakrabarty 
et al. 2001; Murray et al. 2002a). However, one conclusion 
is crystal clear: the X-ray source is not a standard radio 
pulsar (unbeamed or otherwise). 

The basis of our effort is that observationally, all cen- 
tral sources in SNRs known to date, regardless of the band 
of their initial identification (7-ray, X-ray, or radio) ap- 
pear to possess detectable X-ray emission. Theoretically, 
we expect thermal X-ray emission from young neutron 
stars. Thus, on both counts the search for central sources 
in young remnants is very well motivated. However, a 
follow-up program is essential since many other foreground 
sources such as flare stars, young stars, and accreting 
sources and background sources such as AGN dominate 
the source counts (Hertz & Grindlay 1988; Mushotzky 
et al. 2000). Fortunately, the sub-arcsecond spatial res- 
olution of Chandra allow efficient filtering of such contam- 
inating objects. 

To this end, we have identified a sample of SNRs within 
5kpc of the Sun which do not have known radio pulsars 
or PWNe and have not been associated with Type la ex- 
plosions (Table 1. We successfully proposed for a "large" 
Chandra effort in AO-3 to image nine SNRs. This initial 
allocation has been supplemented with additional time in 
AO-4 of Chandra and AO-2 of XMM. Followup of the X- 
ray sources has been undertaken with a plethora of ground 
based telescopes in the optical and near-IR bands (Palo- 
mar 60-inch, Palomar 200-inch, Las Campanas 40-inch, 
ESO 3.5-m, Magellan 6.5-m and Keck 10-m). Here, we re- 
port the first analysis of four SNRs for which the followup 
is now complete. The analysis for the remaining remnants 
will be reported in future papers. 

The organization of the paper is as follows. In § 2 we 
summarize the empirical X-ray properties of the known 
sample of young neutron stars. Such a summary is es- 
sential since the guiding principle of our effort is to place 
our search for central objects against the framework of ex- 
isting classes of sources. Specifically, we are not entirely 
guided by the relatively poorly understood cooling of neu- 
tron stars. Our search has been designed to find objects 
as faint as one tenth of the central X-ray source in Cas- 
siopeia A. In § 3.1 we present a complete sample of cata- 
loged SNRs within 5kpc. Of these, 18 have an identified 
central source or are known to be a composite remnant, 
while three are thought the be the results of Type la SNc. 
The remaining 23 seemingly "hollow" SNRs form our pri- 
mary sample. By "hollow", we refer to SNRs that have 
distinct shells but no obvious indication of central neu- 
tron stars — see Vasisht & Gotthelf (1997). Section 4 has 
a general overview of our observations and analysis tech- 
niques: in § 4.1-4.4 we present a summary of the details of 
our Chandra observations and data reduction, and a like- 
wise global description of the extensive multi-wavelength 
followup is given in § 4.5. We follow this by detailed de- 
scriptions of each of the first four SNRs, its observations, 
and identification of counterparts to its X-ray sources in 
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§§ 5.1, 5.2, 5.3, and 5.4 for SNRs G093. 3+6.9, G315.4-2.3, 
G084.2+0.8, and G127.1+0.5, respectively. It must be ap- 
preciated that complete identification of all X-ray sources 
is essential, given the small sample size. We hope that 
our detailed cataloging will be of help to efforts such as 
ChaMPlane (Grindlay et al. 2003). Finally, in § 6 wc dis- 
cuss what limits our data can place on the existence of 
central sources in the four SNRs, and we conclude in § 7. 

2. X-RAY PROPERTIES OF YOUNG NEUTRON STARS 

The first manifestations of neutron stars associated with 
SNRs were traditional rotation-powered pulsars (such as 
those in the Crab and Vela SNRs). With the advent of 
high-energy missions a number of new classes were dis- 
covered. These include Soft Gamma Repeaters (SGRs; 
for a review see Hurley 2000), Anomalous X-ray Pulsars 
(AXPs; for a review see Mereghetti et al. 2002a), and Com- 
pact Central Objects (CCOs; for a review see Pavlov et al. 
2002). Finally, recent radio surveys have uncovered cen- 
tral Crab-like pulsars with field strengths beyond 10 13 G 
- the so-called high-B pulsars (HBPSRs; Camilo et al. 
2000; Gotthelf et al. 2000). SGRs and AXPs have been 
suggested to be magnetars (Duncan & Thompson 1992; 
Thompson & Duncan 1995): neutron stars with extremely 
strong field strengths, B > 10 15 G. CCO is the generic 
name for a heterogeneous group of X-ray sources emit- 
ting largely unpulsed soft (thermal) emission and lacking 
detectable radio emission. Another possibly related class 
are the Isolated Neutron Stars (INSs; Haberl 2003), also 
called Radio-Quiet Neutron Stars (RQNSs) or Dim Ther- 
mal Neutron Stars (DTNSs). Below we summarize the 
general properties of each of these classes in turn. 

2.1. Radio Pulsars 

Radio pulsars observed in the X-rays often have two- 
component spectra. Here, the thermal component tends 
to be softer (< 0.3 keV) and the power-law component 
harder (r = 1.5-2.5) than those of AXPs; see Becker & 
Trumper (1997), Possenti et al. (2002), and Tables 2 and 
3. A rough relation was originally discovered by Seward & 
Wang (1988) between the X-ray luminosity and the rota- 
tional energy loss rate E: £x,o.i-2.4kcV ~ 10~ 3 -E (Becker 
& Trumper 1997). While this relation has been updated 
for specific classes of neutron stars and different energy 
bands (e.g., Possenti et al. 2002) and has considerable scat- 
ter, it still holds on average . The observed X-ray lumi- 
nosities of radio pulsars then vary between 10 31 ergs s _1 
and 10 37 ergs s _1 , depending on their values of E (and 
through that, their values of P and P). 

While only millisecond pulsars and "old" (> 10 6 yr) pul- 
sars from the sample of Possenti et al. (2002), have values 
of E less than 10 34 ergs s~ 4 , there have been three young 
pulsars discovered recently that are more extreme. These 
arc the HBPSRs PSR J1814-1744 (P = 4.0 s and E = 
4.7 x 10 32 ergs s" 1 ; Camilo et al. 2000), PSR J1847-0130 
(P = 6.7 s and E = 1.7 x 10 32 ergs s _1 ; McLaughlin 
et al. 2003), and PSR J1718-37184 (P = 3.4 s and E = 
1.6 x 10 33 ergs s _1 ; McLaughlin et al. 2003). The first two 
do not have detected X-ray emission (Pivovaroff, Kaspi, & 
Camilo 2000, McLaughlin et al. 2003), while the third does 
have a very faint (L 2 -iokoV ~ 9 x 10 29 ergs s" 1 ) X-ray 



counterpart [the energetic (E > 10 36 ergs s 1 ) HBPSRs 
have brighter X-ray counterparts (Gaensler et al. 2002; 
Gonzalez & Safi-Harb 2003, Helfand, Collins, & Gotthelf 
2003)], consistent with their values of E. 

2.2. Pulsar Wind Nebulae 

Pulsar wind nebulae (PWNc), are bright, centrally con- 
densed nebulae with non-thermal (power-law) X-ray and 
radio spectra often associated with young, energetic pul- 
sars and SNRs (here we refer only to "bubble" PWNe, as 
differentiated by the bow-shock PWNe produced by the 
motion of the pulsars through the ambient medium; for 
reviews, see Chevalier 1998 or Gaensler 2003). The photon 
indices range from 1.3-2.3 (Gotthelf 2003), similar to those 
of pulsars, but they are roughly ~ 10 times as luminous for 
a given E (Gotthelf 2004) and the sizes range from a few 
arcseconds to several arcminutcs. PWNe, both X-ray and 
radio, offer a great advantage over bare radio pulsars for in- 
ferring the existence of neutron stars: they are unbeamed. 
This fact has historically been used in a number of cases 
to infer the existence of energetic pulsars where the pulsar 
had not been seen itself, such as 3C 58 (Becker, Helfand, 
& Szymkowiak 1982, Frail & Moffett 1993), N157B (Wang 
& Gotthelf 1998), and Kes 75 (Becker & Helfand 1984). 

The evolution of a PWN (see van der Swaluw, Downes, 
& Kcegan 2003 for a recent review) begins as it expands su- 
pcrsonically through the shocked ejecta of the SNR. Here, 
the radius of the PWN is oc S 1 / 5 ^/ 5 (assuming that the 
total luminosity of the neutron star is oc E; van der Swaluw 
et al. 2001; van der Swaluw 2003), so it depends only 
weakly on E but more strongly on the remnant age. Even- 
tually, though, the reverse shock of the SNR crashes back 
on the PWN at a time t co \ <~ 10 4 yr (van der Swaluw et al. 
2003, although this depends on the ejecta mass). After this 
occurs , the spherical portion of the PWN continues to ex- 
pand subsonically with its radius oc t 11 / 15 (van der Swaluw 
et al. 2001), while at the same time the pulsar moves away 
from the center of the SNR and begins to form a bow-shock 
nebula (when its motion through the SNR ejecta becomes 
supersonic at ages of > 3000 yr; van der Swaluw et al. 
2003). The passage of the reverse shock should re-heat 
and energize the PWN, causing it to brighten. 

2.3. AXPs & SGRs 

In the X-ray band, traditional AXPs are characterized 
by a two-component spectrum: a power-law with index 
3-4, and a soft blackbody with temperature 0.3-0.7 keV 
(Mereghetti et al. 2002a). While the distances are quite 
uncertain, especially to those not in SNRs, the luminosi- 
ties are relatively constant and > 10 34 ergs s _1 (in the 
1-10 keV band; Mereghetti et al. 2002a), substantially 
greater than their spin-down energy loss rates. These 
sources are characterized by slow spin periods (6-12 s) 
with substantial period derivatives that indicate the pres- 
ence of very strong magnetic fields (B > 10 14 G). 

There are two objects, however, that may indicate 
AXPs can have substantial X-ray variability. First, 
AX J184453-025640 (Vasisht et al. 2000) is a 7 s X- 
ray pulsar which varied by a factor of > 10 in flux but 
whose properties are otherwise found to be fully consis- 
tent with an AXP. More recently, the AXP-candidatc 
XTE J1810-197 was discovered in RXTE data (Ibrahim 
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Table 1 
SNRs Within 5 kpc 



SNR G 


Other Name 


Dist. 

(kpc) 


Distance Method 


Sizc a 
(arcmin) 


Type b 


Central Source/ 
la 


X-ray? d 


Sample 


Refs f 


004.5+6.8 


Keplcr s 


4.5 


Optical expansion/Hi 


3 


S 


Ia?° 






1 


005.4-1.2 


Duck 


5 


H i absorption 


35 


C? 


pulsar? 


yes 




2,3 


006.4-0.1 


W 28 g 


3.0 


OH masers 


42 


S 








4 


nil Q HQ 
Ull .z — U.o 




5 


H i absorption 


A 

4: 


(-1 
\-> 


pulsar 


yes 




C 

O 


013.3-1.3 




2.4 


CO absorption 


70 


s 






XMM-A02 




034.7-0.4 


W44 


2.5 


H i absorption 


35 


c 


pulsar 


yes 




6,7,8 


039.7-2.0 


W50 


5 


VLBI 


120 


? 


SS 433 


yes 




9,10 


053.6-2.2 


3C 400.2 


2.8 


H i association 


33 


s 




CXO-A03 




U04.4 — U.O 


HC 40 


o.o 


CO association 


-1 n 


O 






CYC) AO'} 




065.3+5.7 




0.8 


Optical velocity 


310 


s 






CXO-A04 




069.0+2.7 


CTB 80 


2 


H i absorption 


80 


? 


pulsar 


yes 




11,12 


074.0-8.5 


Cygnus Loop 


0.44 


Optical proper motion 


230 


s 


X-ray source? 




CXO-A04 


13 


078.2+2.1 


7 Cygni 


1.2 


OH/CO association 


60 


s 






XMM-A02 


14 


Uo4.z — U.o 






Hi/CO association 




Q 
O 






OA (J-JWJo 




089.0+4.7 


HB 21 


0.8 


OB association 


120 


s 










093.3+6.9 


DA 530 


3.5 


Hi absorption/X-ray fitting 


27 


s 






CXO-A03 




111.7-2.1 


Cas A 


3.4 


Optical expansion 


5 


s 


CCO 


yes 




15 


114.3+0.3 




3.5 


H i association 


90 


s 


pulsar 


yes 




16,17,18 


lie K i 1 i 
1 ID.O-)- 1 . 1 




A 


H i association 


oU 


Q 
O 










116.9+0.2 


CTB 1 


3.1 


Optical lines 


34 


s 






CXO-A03 




119.5+10.2 


CTA 1 


1.4 


H i association 


90 


s 


7-/X-ray neutron star 






19,20,21 


120.1 + 1.4 


Tycho 


2.4 


Proper motion/shock vel. 


8 


s 


la 






22 


127.1+0.5 


R5 


1.3 


Assoc. with NGC 559 


45 


s 






CXO-A03 




130.7+3.1 


3C 58 


3.2 


H i absorption 


9 


F 


pulsar 


yes 




23,24 


1 1J- 1 ^ 
LOZ. 1 -\-L.O 


HB 3 


2 2 


Interactions with ISJvl 


OU 


c 






VlllM A Pi 9 
AMM -1WJZ 




156.2+5.7 




1.3 


NEI fits 


110 


s 






CXO-A04 




160.8+2.6 


HB9 


< 4 


H i optical vcl. 


140 


s 






CXO-A04 




166.0+4.3 


VRO 42.05.01 


4.5 


H i association 


55 


s 






CXO-A03 




166.2+2.5 




4 


H i interaction 


90 


s 










180 — 1 7 


S147 


1 2 


pulsar dispersion measure 


180 


g 


pulsar 


yes 




oc: ofi 97 


184.6-5.8 


Crab 


2 


proper motion/radial vel. 


7 


F 


pulsar 


yes 




28 


189.1+3.0 


IC 443 


1.5 


Opt. vel. /assoc. with S249 


45 


s 


neutron star 


yes 




29 


203.0+12.0 


Monogcm Ring 


0.3 


O vi absorption/modeling 


1500 


s 


pulsar 


yes 




30 


one c^^n ^ 


MoFiOCCTOS 


1.2 


Optical velocity 


220 


Q 
O 










260.4-3.4 


Puppis A 


2.2 


H i association 


60 


s 


CCO 


yes 




31,32 


263.9-3.3 


Vela 


0.3 


pulsar parallax 


255 


c 


pulsar 


yes 




33 


296.5+10.0 


PKS 1209-51/52 


2.1 


H i association 


90 


s 


CCO/X-ray pulsar 


yes 




34,35 


309.8+0.0 




3.6 


H i absorption 


25 


s 






CXO-A03 




315.4-2.3 


RCW 86 


2.8 


Optical lines 


42 


s 






CXO-A03 


36 


327.4+0.4 


Kes 27 


4.3 


Hi absorption/interact. 


30 


s 










327.6+14.6 


SN 1006 


2.2 


Spectra/proper motion 


30 


s 


la 






37,38 


330.0+15.0 


Lupus Loop 


0.8 


Hi 


180 


s 






CXO-AOi 




332.4-0.4 


RCW 103 


3.3 


H i absorption 


10 


s 


CCO 


yes 




39 


343.1-2.3 




2.5 


H i absorption 


32 


c 


pulsar 


yes 




40,41 


354.1+0.1 




5 


recombination lines 


15 


c 


pulsar? 


no 




42,43 



a Major axis. 

b Typcs are: S (shell), C (composite), and F (filled, or PWN), as determined by Green (2001). A question mark indicates that the type 
is poorly determined. 

c Indicatcs if a central source is known, or if the SNR is thought to be of type la. If a central source is known and it falls into one of the 
classes from § 2, it is labeled accordingly. 

d Indicates whether the central source has been detected in X-rays. See Table 2. 
G Rcfers to X-ray samples of shell SNRs described in § 4.1. 

f References deal only with the central source or la classification. General SNR properties were taken from Green (2001) and references 
therein. 

g Has already been observed with Chandra. 
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34: Hclfand & Becker (1984); 35: Zavlin et al. (2000); 36: Vink et al. (2000); 37: Fcscn ct al. (1988); 38: Allen ct al. (2001); 39: Gotthelf, 
Petre, & Hwang (1997); 40: McAdam, Osborne, & Parkinson (1993); 41: Becker, Brazier, & Triimper (1995); 42: Frail ct al. (1994); 43: 
Becker & Triimper (1997). 
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Table 2 



X-ray Properties of Central Sources from Table 1 



Source 


SNR G 




r 


log 10 L x a 


X-ray 


Age b 


Refs 






(keV) 




(ergs s -1 ) 


PWN? 


(kyr) 




PSR B1757-24 


005.4-1.2 




1.6 


33.0 


Ycs c 


> 39 


1,2 


PSR J1811-1925 


011.2-0.3 




1.7 


33.7 


Yes 


1.6 


3,4,5 


PSR B1853+01 


034.7-0.4 




1.3 


31.2 


Yes 


20 


6 


SS 433 c 


039.7-2.0 




0.7 


35 




5-40 


7 


PSR B1951+32 


069.0+2.7 




1.6 


32.9 


Yes 


64 


8 


CXO J232327.8+584842 


111.7-2.1 


0.7 


3.0 


33.5/34.7 h 


No 


0.3 


9 


PSR B2334+61 


114.3+0.3 




2.0 d 


31.7 


... g 


41 


10 


RX J0007.0+7302 


119.5+10.2 


0.14 


1.5 


31.2 


Yes 


13 


11 


PSR J0205+6449 


130.7+3.1 




1.7 


32.2 


Yes 


0.8 


12,13 


PSR J0538+2817 


180.0-1.7 


0.16 




32.9 


Yes 


30 


14,15 


PSR B0531+21 


184.6-5.8 




1.6 


35.8 


Yes 


1.0 


16 


CXO J061705.3+222127 


189.1+3.0 


0.7 




31.3 


Yes 


30 


17 


PSR B0656+14 


203.0+12.0 


0.07+0.14 


1.5 


31.1 


No 


100 


18,19,20 


RX J0822-4300 


260.4-3.4 


0.4 




33.5 


No 


3.7 


21,22 


PSR B0833-45 


263.9-3.3 




2.5 


32.5 


Yes 


11 


23 


IE 1207.4-5209 


296.5+10.0 


0.26 




33.1 


No 


7 


24 


IE 161348-5055° 


332.4-0.4 


0.6 




32.1-33.9 


No 


2 


25 


PSR B1706-44 


343.1-2.3 


0.14 


2.0 


32.6 


Yes° 


18 


26 


PSR B1727-33 


354.1+0.1 






< 32.6/< 32.5 


... s 


15 


16 



References.— 1: Kaspi ct al. (2001); 2: Gaensler & Frail (2000); 3: Torii et al. (1997); 4: Roberts et al. 
(2003); 5: Stephenson & Green (2002); 6: Pctrc ct al. (2002); 7: Kotani et al. (1996); 8: Migliazzo ct al. (2002); 
9: Mereghetti, Ticngo, & Israel (2002b); 10: Becker ct al. (1996); 11: Slanc ct al. (2004); 12: Slane et al. 
(2002); 13: Stephenson (1971); 14: Romani & Ng (2003); 15: Kramer ct al. (2003); 16: Becker & Triimper 
(1997); 17: Olbert et al. (2001); 18: Marshall & Schulz (2002); 19: Briskcn ct al. (2003); 20: Greiveldinger 
ct al. (1996); 21: Zavlin et al. (1999); 22: Pavlov et al. (2002); 23: Pavlov ct al. (2001b); 24: Sanwal ct al. 
(2002); 25: Gotthelf, Petre, & Vasisht (1999); 26: Gotthelf, Halpern, & Dodson (2002). In addition, general 
pulsar data have been taken from Hobbs & Manchester (2003). 

a Luminosity for only the point-source in the 0.5-2.0 keV band, assuming the distance from Tabic 1. Upper 
limits to the luminosity are given for a blackbody with fcToo = 0.25 keV and for a power-law with T = 2.2. 

b The best estimate of the age of the SNR if known, otherwise the pulsar spin-down age P/2P. 

c Possibly not an isolated neutron star. 

d Assumed. 

e The X-ray PWNe here are significantly fainter compared to the pulsars than for other sources (Kaspi et al. 
2001; Gotthelf et al. 2002). 

g The current X-ray data do not sufficiently constrain the existence of a nebula. 

h The current X-ray data do not constrain the spectrum: either a blackbody or power-law model is possible. 
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Table 3 

Properties of Rotation-Powered Pulsars Associated with SNRs 



PSR 


SNR 


P 


T a 


login E 


D 


logio L x b 


Refs. 






(ms) 


(kyr) 


(ergs s 1 ) 


(kpc) 


/ — 1 \ 
(ergs s ) 




10205+6449 


ci 30 7+3 1 


66 


5.4 


37.4 


3.2 


32.2 




J0537-6910 


N157B C 


16 


5.0 


38.7 


49.4 


38.3 


1,2 


J0538+2817 


G180.0-1.7 


143 


620 


34.7 


1.8 


32.9 




B0531+21 


G184.6-5.8 


33 


1.2 


38.7 


2.0 


35.8 




B0540-69 


0540-69. 3 C 


50 


1.7 


38.2 


49.4 


36.3 


2,3 


B0833— 45 


G263 9—3 3 


89 


11 


36.8 


0.3 


32.5 




J1016-5857 


G284.3-1.8 


107 


21 


36.5 


3 


32.5 


4 


J1119-6127 


G292.2-0.5 


408 


1.6 


36.4 


6 


32.5 


5 


J1124-5916 


G292.0+1.8 


135 


2.9 


37.1 


> 6 


33.0 


6,7 


B1338-62 


G308. 8-0.1 


193 


12 


36.1 


9 


32.8 


8 


B1509-58 


G320.4-1.2 


151 


1.6 


37.3 


5.2 


34.1 


9,10 


B1643-43 


G341. 2+0.9 


232 


33 


35.5 


7 






B1706-44 


G343. 1-2.3 


102 


18 


36.5 


2.5 


32.6 




B1727-33 


G354.1+0.1 


139 


26 


36.1 


5 


<32.6 




B1757-24 


G005.4-1.2 


125 


16 


36.4 


5 


33.0 




J1811-1925 


G011.2-0.3 


65 


23 


36.8 


5 


33.7 




J1846-0258 


G029.7-0.3 


324 


0.72 


36.9 


19 


35.4 


11 


B1853+01 


G034.7-0.4 


267 


20 


35.6 


2.5 


31.2 




J1930+1852 


G054.1+0.3 


136 


2.9 


37.1 


5 


33.3 


12 


B1951+32 


G069.0+2.7 


40 


110 


36.6 


2 


35 




J2229+6114 


G106.6+2.9 


52 


11 


37.3 


3 


32.8 


13,14 


B2334+61 


G114.3+0.3 


495 


41 


34.8 


3.5 


31.7 





Note. — Pulsar-SNR associations are largely drawn from Kaspi & Helfand (2002). 
Characteristic age r = P/2P. 

b X-ray luminosity in the 0.5-2.0 keV band. Upper limits to the luminosity are given for 
a power-law with r = 2.2. 

c In the Large Magellanic Cloud. 

References.— 1: Wang et al. (2001); 2: Gotthelf (2003); 3: Gotthelf & Wang (2000); 4: 
Camilo et al. (2001); 5: Gonzalez & Safi-Harb (2003); 6: Hughes et al. (2001); 7: Gaensler 
& Wallace (2003); 8: Gaensler, Kaspi, & Manchester (2003); 9: Greiveldinger et al. (1995); 
10: Gaensler et al. (2002); 11: Helfand ct al. (2003); 12: Camilo et al. (2002a); 13: Kothes, 
Uyaniker, & Pineault (2001); 14: Halpcrn ct al. (2001) . Also see references from Tabs. 1 
and 2. 
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et al. 2003); Chandra and XMM data confirm the pulsa- 
tions present in the RXTE data and allow for comparison 
with archival Einstein, ROSAT, and ASCA data where it 
is a factor of s» 100 fainter (Gotthelf et al. 2003). In the 
bright (current) state, XTE J1810— 197 has an absorbed 
X-ray flux (0.5-10 keV) of ps 4 x lO" 11 ergs citT 2 s~\ 
while ROSAT data from 1993 has a flux of w 5 x 
10 -13 ergs cm~ 2 s _1 . Converting the fluxes to luminosi- 
ties is uncertain due to the largely unconstrained dis- 
tance, but Gotthelf et al. (2003) determine an upper limit 
of 5 kpc. With this, the current luminosity in the 0.5- 
2 keV range is < 3 x 10 35 ergs s _1 , similar to other AXPs, 
but the "quiescent" luminosity is ~ 10 33 ergs s _1 . While 
less certain, AX J184453— 025640 shows a roughly similar 
range of luminosities. 

SGRs have roughly similar quiescent luminosities (Hur- 
ley 2000), although their spectra are somewhat harder. 
However, none of the SGRs is firmly associated with a 
SNR (Gaensler et al. 2001). There is now mounting ob- 
servational evidence that SGRs and AXPs are related ob- 
jects (Gavriil, Kaspi, & Woods 2002, Kulkarni et al. 2003; 
Kaspi et al. 2003), confirming the hypothesis of Thomp- 
son & Duncan (1996). In what follows, we do not treat 
the SGRs as separate objects since their quiescent X-ray 
properties are sufficiently similar to those of AXPs. 

2.4. Cooling Radio-Quiet Neutron Stars/Compact 
Central Objects 

Regardless of emission at other wavelengths, young neu- 
tron stars should have thermal X-ray emission from their 
surfaces as they cool. Indeed, it is likely that ROSAT 
has detected a number of nearby, ~ 10 6 yr (too old to be 
in SNRs) neutron stars (Treves et al. 2000; Haberl 2003) 
with no detectable radio emission (Kaplan, Kulkarni, & 
van Kerkwijk 2003a, Kaplan et al. 2003b; Johnston 2003) 
— hence the name RQNS. The ROSAT sources have tem- 
peratures of < 100 eV. This thermal emission is almost 
certain to be present at some level, regardless of whatever 
other processes are occurring (radio emission, accretion, 
etc.), but the exact level of emission depends on the mass 
of the neutron star and on the presence or absence of exotic 
particles (pions, kaons, hyperons, free quarks, etc.) and/or 
processes (i.e. direct Urea cooling) in the core (Kaminker, 
Yakovlev, & Gnedin 2002, Tsuruta et al. 2002). 

A number of superficially similar sources known as 
CCOs have been discovered in SNRs, such as the sources 
in Cas A and Puppis A. The temperatures of CCOs should 
be higher than those of the field RQNSs by a factor of 5 
or so, depending on cooling physics (e.g., Page 1998). For 
some of these sources, it is likely that there is little if any 
radio emission due to the low values of E (Kaplan et al. 
2002a), but for others the lack of such emission may just be 
beaming effect (i.e., as with Geminga), implying beaming 
fractions of ~ 50% (Brazier & Johnston 1999). The radii 
of the blackbody fits to the CCOs in Table 2 are typically 
less than the 10-15 km expected for a neutron star (pre- 
sumably since a blackbody, while providing an adequate 
spectral fit, does not actually represent the emission of the 
surface) — values of R^ « 1-3 km are common. 

3. SURVEY DESIGN 

3.1. A Volume- Limited Sample of Shell SNRs 



The success of this effort hinges upon defining an objec- 
tively constructed sample so that strong conclusions can 
be drawn not only from detections but also non-detections. 
Bearing this in mind, we identified all of the SNRs 6 that 
are at a distance of less than 5 kpc as determined from a 
reliable distance determination (e.g. derived from the kine- 
matic velocity of associated line emission or a parallax, 
rather than from Nn or the £ — D relation) . This sample 
is comprised of the 45 SNRs listed in Table 1. This is not 
an entirely complete sample, as illustrated by the void in 
the third quadrant in Figure 1 and the relative paucity of 
distant SNRs toward the Galactic center, but the criteria 
for inclusion in the sample (detected SNR with a robust 
distance) should not be correlated with the properties of 
the central objects. We estimate that we cover ~ 15% 
of the Galactic molecular gas, and hence sites of massive 
star formation, with our distance criteria (based on Dame 
1993). For the X-ray observations we eliminated all SNRs 
that were of Type la or those that are already associated 
with a NS and/or central synchrotron nebula 7 — i.e. we 
only include those that are type S in Table 1. 

The shell remnants are the major remaining sample 
where the neutron star population has not been systemat- 
ically assessed. For the SNRs with PWNe, some estimate 
as to the neutron star's properties can be made even if 
the source itself has not been detected, but the shell rem- 
nants permit no such estimation. Therefore, these shell- 

6 Drawn from the Galactic SNR Catalog (Green 2000). Since we 
constructed the table, the catalog has been updated (Green 2001), 
and it is this list and more recent references that we use to determine 
the properties of the SNRs. 

7 We assume that sources like IC 443, which has a synchrotron nebula 
but where pulsations have not yet been detected, still do have central 
neutron stars. 



180 




O^GC 



270 

Fig. 1. — Galactic distribution of SNRs from Table 1. Galactic 
longitude I (degrees) is plotted against distance (kpc), with the Sun 
at the origin and the Galactic center to the right. Black squares are 
the la SNRs, blue diamonds are associated with radio pulsars, green 
asterisks are associated with other types of neutron stars (RQNS, 
etc.), and the red circles are the hollow SNRs. The four filled circles 
are the SNRs from this paper with detailed analyses. 
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type SNRs are the subject of our survey. Below we first the 
status of neutron star velocities, an issue that affect our 
survey design. We then continue to describe the survey 
itself. 

3.2. Neutron Star Velocities 

Radio pulsars have high space velocities, among the 
highest in the Galaxy, as measured from timing, scintil- 
lation, and interferometry. Various authors (e.g., Hansen 
& Phinney 1997; Lorimcr et al. 1997; Cordes & Chernoff 
1998; Arzoumanian et al. 2002) have modeled the distribu- 
tion slightly differently, but they all agree that a substan- 
tial number of pulsars move with speeds > 300 km s _1 , 
while i=s 90% have speeds < 700 km s" 1 . 

The space velocities of neutron stars can also be inferred 
from their offsets from the centers of associated supernova 
remnants. Such an approach demands reliable estimates 
for the distance and age of the system under considera- 
tion. It also assumes that the geometric center of a SNR 
is easily identifiable, and that this center corresponds to 
the site of the supernova. With these caveats in mind, such 
an analysis potentially provides a direct measurement of 
the neutron star velocity distribution, free from the many 
selection effects associated with proper motion studies. 

Frail, Goss, & Whiteoak (1994) carried out a detailed 
study of the offsets of radio pulsars from the centers 
of SNRs. From a sample of 15 pulsar/SNR associa- 
tions, Frail et al. (1994) concluded that young radio 
pulsars are born with projected space velocities ranging 
between 50 km s _1 and 1000 km s" 1 , with a mean of 
~500 km s" 1 . While at that time this distribution repre- 
sented a somewhat more rapidly-moving population than 
that which had been inferred from proper motion studies, 
it agrees well with the more recent determinations dis- 
cussed above. 

However, the preceding discussion assumes that pulsars 
were born near the centers of what are now SNRs. Some 
authors dispute that this is always the case (Gvaramadze 
2002), suggesting that if the SN progenitors have large 
space velocities they could evacuate a bubble with their 
winds, move toward the edges of those bubbles, and then 
explode (in a so-called off-center cavity explosion). This 
would make the origins of neutron stars closer to the edges 
of SNRs. While an interesting possibility for a few sources, 
the large number of associations where the neutron star is 
close to the center of the SNR mean that this hypothesis 
cannot work for the majority of sources. 

It is not yet known whether other populations of neu- 
tron star have different velocity distributions from that 
seen for radio pulsars. The location of several SGRs 
on or beyond the perimeters of SNRs originally led to 
the suggestion that SGRs had very high space veloci- 
ties, v± <~ 1000 — 2000 km s~\ as might result from 
the anisotropic neutrino emission associated with the for- 
mation of such highly-magnetized objects (e.g., Duncan 
& Thompson 1992). However, it has since been argued 
that many of these SGR/SNR associations are spurious 
(Lorimer & Xilouris 2000; Gaensler et al. 2001), in which 
case these inferred velocities are not valid. On the other 
hand, several AXPs have convincing associations with 
young SNRs. In all these cases, the AXP lies close to 
the geometric center of its SNR, implying projected space 
velocities for this population < 500 km s~ (Gaensler, 



Gotthelf, & Vasisht 1999, Gaensler et al. 2001), consistent 
with the velocities of radio pulsars. The emerging and still 
enigmatic class of central compact objects (CCOs) are also 
centrally located in young SNRs (see Pavlov et al. 2002 for 
a review), and most likely also have relatively low space 
velocities. 

What little is known about the velocities of older neu- 
tron stars that are not radio pulsars roughly agrees with 
the situation for radio pulsars: the velocities are high, 
> 100 km s -1 . Specifically, the velocities of the INSs 
RX J1856.5-3754 and RX J0720.4-3125 are both w 
200 km s _1 (Kaplan, van Kerkwijk, & Anderson 2002b, 
Motch, Zavlin, & Haberl 2003). 

4. OBSERVATIONS AND DATA ANALYSIS 

In this section we give an overview of the analysis pro- 
cedure that we used for the different SNRs. We start by 
describing the splitting of the SNRs into observationally- 
based sub-samples, of which the Chandra ACIS sub- 
sample is the major component discussed here (§ 4.1). We 
then describe the analysis of the Chandra data that were 
used to identify potential compact objects (§ 4.2). Fi- 
nally, we describe the motivation for and basic analysis of 
the optical and infrared followup observations that were 
used to reject contaminating foreground and background 
X-ray sources (§ 4.5). Following this section we present 
the actual detailed analyses of the four SNRs in this paper 
(§ 5). 

4.1. Chandra ACIS Survey 

We defined three sub-samples among the 23 SNRs that 
had no central sources from Table 1 . The sub-samples were 
defined largely by size so that the X-ray observations have 
a good chance of identifying the central compact source. 
Our primary sub-sample was designed for the Chandra 
ACIS-I detector, with its 16' ficld-of-view. We parame- 
terize as follows: a neutron star has a typical transverse 
velocity of lOOuioo km s _1 , distance d in kpc, and an age 
of 10 4 i 4 yr. To ensure that a NS lies within 8' of its SNRs 
geometric center (and so will fall on the ACIS-I array; see 
Figs. 2-5 for illustrations of the ACIS field-of-view), we 
require t 4 < 2.27d/wioo- For a Sedov-phase remnant, we 
have U = 7 x lO _3 0<ff 7 ~ 1/2 , where 9 arcmin is the SNR di- 
ameter and 10 7 T7 K is the shock temperature. Our limit 
is thus 9 < 324Tj /2 v w \ ) . We expect T 7 w 1 for a broad 
range of SNR ages; given the somewhat weak dependence 
on T we adopt T7 = 1 in this calculation. A conserva- 
tive constraint on pulsar velocities is «ioo < 7. This then 
yields the condition 9 < 46 arcmin. There arc 14 SNRs 
that meet this criterion, of which 3 (Cas A, IC 443, and 
Kepler) have already been observed by Chandra. This then 
leaves 11 SNRs for further Chandra observations, although 
SNR G327.4+0.4 only recently had a distance determina- 
tion and was not included in the original Chandra sample. 
We therefore proposed for ACIS-I observations of the 10 
remaining SNRs (identified by CXO- A03 in Tab. 1), and 
were awarded observations of 9 (the tenth, G006.4— 0.1, 
was awarded to J. Rho in another AO-3 proposal). 

We then constructed a sample of the larger SNRs for 
XMM, and one of the largest SNRs using the literature 
and the ROSAT All-Sky Survey Bright Source Catalog 
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Fig. 3.— MOST (Whiteoak & Green 1996) 833 MHz radio image 
of SNR G315.4-2.3, showing the placement of the ACIS-I detector. FlG - 5 -~ DRAO 1.4 GHz radio image of SNR G127.1+0.5, show- 

The bright region toward the south-west is the nebula RCW 86 ln § the placement of the ACIS-I detector. 

(Rodgcrs, Campbell, & Whiteoak 1960). The contours are from 
ROSAT PSPC data (sequence RP500078). 
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(BSC) (Voges et al. 1996), coupled with Chandra snap- 
shot images. These samples are identified by XMM-A02 
and CXO-AOA in Table 1 (our XMM proposal was for 
seven SNRs, of which we were only allocated time for the 
three that are identified as "XMM-A02'). We defer de- 
tailed discussion of these samples to later papers. 

To determine exposure times for these sources, we ex- 
amined the types of neutron stars found in the SNRs from 
Table 1 (see also Chakrabarty et al. 2001 and § 2). These 
neutron stars are listed in Table 2. We see two groups 
among them: those with non-thermal spectra (charac- 
terized by a power-law with photon index T ~ 1.6) and 
those with thermal spectra (characterized by a blackbody 
with temperature fcToo ~ 0.5 keV). Among the shell- type 
SNRs in Table 1, the thermal sources predominate, and 
these also produce the lowest X-ray count-rates. There- 
fore we computed exposure times for a thermal source with 
fcToo = 0.25 keV (toward the low-end of those in Tab. 2) 
and a bolometric luminosity Lboi = 10 32 ergs s _1 , a factor 
of ~ 10 lower than those of most of the thermal sources. 
For the column densities JVh we used the best available es- 
timates from the literature. For these shell SNRs, there is 
very little contribution from the SNR itself in the interior, 
especially with the resolution of Chandra, so this was not 
an issue, although we calculate the expected background 
contribution, using surface-brightnesses compiled from the 
literature. The exposure times are those that should de- 
tect a source with a prototypical spectrum but a factor 
of 10 less luminous than normal with a significance of at 
least 5 a above the background. 

4.2. X-ray Point-Source Detection 

To analyze the data, we first reprocessed the raw (level 
1) event data to take advantage of updates in the Chandra 
calibration since the data were first taken. Specifically, the 
reprocessing included a correction for charge transfer inef- 
ficiency (CTI; Townsley et al. 2002) 8 , and we removed the 
±0.5 pixel randomization added to the events. We did not 
include any correction for the degradation of the quantum 
efficiency (QE) of the ACIS detectors 9 , as at no point did 
we do a complete spectral analysis that would have used 
the available correction techniques 10 . These corrections 
are minor (« 10% year -1 ) and are beyond our level of ac- 
curacy. We selected only events that have the "standard" 
ASCA grades (0, 2, 3, 4, and 6). After generating a new 
level 2 event file, we corrected the data for Chandra aspect 
errors 11 : for example, the change for SNR G093.3+6.9 
was — 0'.'13 in Right Ascension and 0'.'32 in Declination. 
Smoothed images of the data are shown in Figure 6-9. 

Starting from this level 2 event file, we processed the 
data and extracted sources much as was done in the Chan- 
dra Deep Field (CDF; Brandt et al. 2001) and the 82-ksec 
observation of the Orion Nebula (Feigelson et al. 2002). 
First we separated the level 2 events by energy into three 
bands: 0.3-2.0 keV (L), 2.0-8.0 keV (H), and 0.3-8.0 keV 
(A). These bands are slightly modified from those of the 
CDF and Orion, as our potential sources are softer. We 
then separated the data by CCD, using only the four 

8 Following http://asc.harvard.edu/ciao/threads/acisapplycti/. 

9 See http://cxc.harvard.edu/cal/Acis/Cal_prods/qeDeg/. 

10 See http://cxc.harvard.edu/ciao/threads/apply_acisabs/. 

11 Following http://asc.harvard.edu/cal/ASPECT/fix_ofTset/fix_ofTset.cgi. 




Fig. 6.— Smoothed ACIS-I image (0.3-8.0 keV) oT 
SNR G093. 3+6.9. North is up, and East is to the left. The bright- 
ness is scaled proportional to the logarithm of the counts in 2" bins, 
smoothed with a Gaussian filter. The sources from Table 4 are 
labeled. 




Fig. 7.— Smoothed ACIS-I image (0.3-8.0 keV) oT 
SNR G315.4-2.3. North is up, and East is to the left. The bright- 
ness is scaled proportional to the logarithm of the counts in 2" bins, 
smoothed with a Gaussian filter. The sources from Table 5 are 
labeled. The diffuse emission to the West and South- West is the 
RCW 86 complex. 
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Fig. 8.— Smoothed ACIS-I image (0.3-8.0 keV) of 
SNR G084.2-0.8. North is up, and East is to the left. The bright- 
ness is scaled proportional to the logarithm of the counts in 2" bins, 
smoothed with a Gaussian filter. The sources from Table 6 are 
labeled. 




Fig. 9.— Smoothed ACIS-I image (0.3-8.0 keV) of 
SNR G127.1+0.5. North is up, and East is to the left. The bright- 
ness is scaled proportional to the logarithm of the counts in 2" bins, 
smoothed with a Gaussian filter. The sources from Table 7 are 
labeled. 
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ACIS-I CCDs. This then gave twelve separate event files 
(3 bands x 4 CCDs). 

For each of these event files, we generated instrument 
and exposure maps, following standard CIAO threads 12 . 
Creating an exposure map requires a source model: for 
this model, we assumed a blackbody with a temperature 
of kToo — 0.25 keV (similar to that of the Cas A central 
point source) and a column density appropriate for each 
SNR (§§ 5.1-5.4). 

Then, for each of the event files (now with exposure 
maps), we ran wavdetect (Freeman et al. 2002), using 
wavelet scales of 1, \/2, 2, \/8, . . • , 16 pixels and a nominal 
energy of 1.5 keV. The significance threshold was 10~ 6 . 

We then merged the wavdetect source files for further 
analysis. First we merged by CCD, creating a source 
list for each band. Then, to merge the data from dif- 
ferent bands (which are spatially coincident), we consid- 
ered sources identical if the positions (as determined by 
wavdetect) differed by less than 2"5, for sources at off-axis 
angles of < 6', or differed by less than 4", for sources at off- 
axis angles of > 6'. While these are large tolerances given 
the typical astrometric uncertainty of wavdetect (0'.'6 for 
on-axis sources, approaching 1'.'5 for off-axis sources), the 
source density is so low (~ 0.2 arcmin~ 2 for this observa- 
tion) that the number of false matches is negligible. Even 
so, only 5 sources had match-radii of > 1'.'5, and most 
were < 0'.'5. After manually examining the list of sources, 
we then removed those that appeared spurious. We then 
limited the source list to those sources with > 10 counts. 
This limit allows us to have enough photons that the po- 
sition is well determined (critical when trying to identify 
counterparts) and that an estimate can be made of the 
hardness ratio. Given that our exposure times were cal- 
culated such that a low-luminosity neutron star would be 
detected with > 25 counts, the 10-count limit is conserva- 
tive. The final merged source list contains 12-18 sources, 
depending on the SNR. The sources are all consistent with 
being point-sources. 

For each source, we then performed additional aperture 
photometry. This allowed us to (1) use knowledge of the 
psf size in determination of the source counts and (2) ob- 
tain source counts for sources that were detected in only 
one or two bands. The radii of the photometric apertures 
were determined from the analytic fit to the 95% encir- 
cled energy radius given in Feigelson et al. (2002). We 
measured the number of counts in these apertures for each 
band and subtracted the number of counts in a background 
annulus extending from 2-3 times the 95% encircled en- 
ergy radius to determine the net number of counts. To 
aid in comparison with the CDF/Orion data sets, we also 
extracted photons in the more standard 0.5-2.0 keV band. 
The final source data are presented in Tables 4-7. We plot 
the L-band counts vs. the H-band counts for the detected 
sources in Figures 10-13. 

4.2.1. Nomenclature 

Herein, for convenience, we label the CTicmdra-detected 
X-ray sources by their field identification — for exam- 
ple, SNR G093.3+6.9:5 refers to the 5th source in the 
SNR G093.3+6.9 field. This is not meant to imply that all 
of these sources are associated with the SNRs — the vast 



http: / /asc. harvard.edu/ciao/threads/expmap_acis_single/. 
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Table 4 



X-ray Sources in SNR G093.3+6.9 



ID a 


CXOU J b 


a 


S 


rgo" 


AR d 


CountSL 


Countsn 


CountSA 


Counts 


HR L ,H° 






(J2000) 


(") 


(') 








(0.5 - 2.0 kcV) 




1 


205222.8+552343 


on c o ■ i ■ i on 

20:52:22.89 


+55:23:43.7 


0.5 


3.5 


310(20) 


47f71 


360(20) 


290(20) 


— 0.73(0.04) 


2 


205230 9+551437 


on co on no 

20:52:30.98 


+55: 14:37.5 


0.6 


6.4 


270f 20 s ! 


127(12) 


390(20) 




— 35f0 05*1 

yj . o { yj . y.f o \ 


4 


205225.8+552741 


20:52:25.82 


+55:27:41.6 


1.3 


7.4 


47(7) 


48(7) 


93.9(10.2) 


44(7) 


0.01(0.11) 


5 


205222.0+551516 


20:52:22.01 


+55:15:16.9 


0.8 


5.3 


20(5) 


13(4) 


33(6) 


19(4) 


-0.20(0.18) 




2052S1 '34-552n'31 


20:51:31.37 


+55:20:31.6 


1.4 


6.1 






32(6) 


1 3f41 


n i i q\ 


7 


205250.2+551606 


20:52:50.20 


+55:16:06.6 


1.5 


6.8 


23(5) 


2.2(1.8) 


25(5) 


22(5) 


-0.83(0.14) 


8 


205242.1+551409 


20:52:42.19 


+55:14:09.8 


1.6 


7.5 


17(4) 


8.0(3.2) 


25(5) 


16(4) 


-0.36(0.21) 


9 


205226.4+551746 


20:52:36.48 


+55:17:46.0 


0.7 


4.2 


11(3) 


7.1(2.8) 


18(4) 


11(3) 


-0.19(0.24) 


10 


205139.8+552553 


20:51:39.85 


+55:25:53.4 


1.7 


7.3 


11(4) 


11(3) 


22(5) 


9.3(3.5) 


-0.02(0.23) 


13 


205242.2+552607 


20:52:42.26 


+55:26:07.1 


1.8 


6.9 


9.3(3.3) 


13(4) 


22(5) 


9.5(3.3) 


0.16(0.23) 


14 


205205.9+551758 


20:52:05.96 


+55:17:58.5 


0.7 


2.8 


2.1(1.5) 


9.2(3.2) 


11(3) 


2.1(1.5) 


0.63(0.24) 


15 


205209.8+551821 


20:52:09.87 


+55:18:21.2 


0.8 


2.2 


7.4(2.8) 


4.0(2.1) 


11(3) 


7.4(2.8) 


-0.30(0.30) 


16 


205250.0+552025 


20:52:50.09 


+55:20:25.6 


1.1 


5.1 


10(3) 


1.4(1.5) 


12(4) 


10(3) 


-0.77(0.24) 


17 


205212.6+551854 


20:52:12.65 


+55:18:54.4 


0.6 


1.6 


3.2(1.8) 


6.3(2.6) 


10(3) 


3.2(1.8) 


0.33(0.31) 


19 


205236.8+551528 


20:52:36.89 


+55:15:28.8 


1.4 


6.0 


4.3(2.4) 


8.4(3.0) 


12(4) 


3.2(2.1) 


0.33(0.29) 


24 


205152.2+552602 


20:51:52.23 


+55:26:02.7 


1.8 


6.4 


7.7(3.0) 


4.5(2.4) 


12(4) 


7.1(2.8) 


-0.26(0.30) 


26 


205248.3+551422 


20:52:48.32 


+55:14:22.3 


3.1 


7.8 


5.9(2.8) 


5.3(2.6) 


12(4) 


5.9(2.8) 


-0.05(0.34) 


30 


205129.7+551548 


20:51:29.72 


+55:15:48.5 


3.0 


7.9 


5.1(2.6) 


9.0(3.3) 


14(4) 


4.0(2.4) 


0.28(0.29) 



internal identifier of the form SNR G093.3+6.9: N. 

b Official IAU name. 

c Approximate 90% confidence radius. 

d Angle from the center of the SNR. 

c Hardness ratio, computed according to (Ch — C'l)/(C , h + Cl), where C is the number of counts in a band. 



Note. — Quantities in parentheses are l-<7 uncertainties. Here and in Tables 5-7, Counts^ — Counts^ + Counts^f , but the three columns have 
been rounded separately. 



Table 5 

X-ray Sources in SNR G315.4-2.3 



ID a 


CXOU J b 


a 


S 




AR d 


CountSL 


CountSH 


CountSA 


Counts 


HRl,h c 






(J2000) 


(") 


(') 








(0.5 - 2.0 koV) 


1 


144319.3-622804 


14:43:19.31 


-62:28:04.2 


0.7 


3.4 


20(5) 


25(5) 


45(7) 


19(4) 


0.11(0.15) 


2 


144333.7-622928 


14:43:33.74 


-62:29:28.2 


1.1 


5.2 


20(5) 


1.6(1.5) 


22(5) 


20(5) 


-0.85(0.13) 


3 


144151.5-622833 


14:41:51.59 


-62:28:33.5 


2.0 


6.8 


11(4) 


10(3) 


22(5) 


12(4) 


-0.05(0.23) 


5 


144219.5-622834 


14:42:19.54 


-62:28:34.5 


1.1 


3.5 


2.9(1.8) 


8.4(3) 


11(3) 


2.9(1.8) 


0.49(0.28) 


6 


144320.6-623308 


14:43:20.63 


-62:33:08.6 


2.0 


6.0 


12(4) 


-0.5(0.2) 


11(4) 


12(4) 


-1.09(0.04) 


8 


144341.0-623138 


14:43:41.05 


-62:31:38.3 


1.9 


6.8 


13(4) 


2.4(1.8) 


16(4) 


14(4) 


-0.69(0.21) 


9 


144346.3-622413 


14:43:46.37 


-62:24:13.1 


3.4 


7.7 


16(5) 


2.2(1.8) 


16(5) 


14(4) 


-0.77(0.19) 


12 


144315.3-622128 


14:43:15.36 


-62:21:28.8 


1.6 


7.5 


8.0(3.2) 


6.1(2.8) 


14(4) 


7.2(3) 


-0.13(0.30) 


13 


144350.4-623040 


14:43:50.45 


-62:30:40.0 


2.2 


7.4 


12(4) 


3.2(2.1) 


15(4) 


12(4) 


-0.57(0.25) 


17 


144322.1-623219 


14:43:22.12 


-62:32:19.3 


1.7 


5.4 


4.3(2.4) 


10(3) 


15(4) 


4.3(2.4) 


0.41(0.27) 


19 


144210.9-622202 


14:42:10.94 


-62:22:02.8 


2.7 


7.7 


10(4) 


3.5(2.1) 


14(4) 


8.8(3.3) 


-0.50(0.26) 


21 


144252.1-622107 


14:42:52.16 


-62:21:07.1 


2.4 


7.2 


10(4) 


5.8(2.6) 


17(4) 


9.8(3.5) 


-0.28(0.26) 


23 


144213.1-623220 


14:42:13.16 


-62:32:20.4 


1.5 


5.8 


9.2(3.2) 


3.7(2.1) 


12(4) 


9.0(3.2) 


-0.43(0.27) 


26 


144207.1-622204 


14:42:07.16 


-62:22:04.3 


3.8 


8.0 


7.5(3.2) 


4.0(2.4) 


12(4) 


8.0(3.2) 


-0.30(0.33) 


28 


144320.0-622138 


14:43:20.06 


-62:21:38.7 


1.9 


7.5 


11(4) 


0.6(1.5) 


12(4) 


11(4) 


-0.90(0.24) 



a Internal identifier of the form SNR G315.4-2.3: TV. 

b Official IAU name. 

c Approximatc 90% confidence radius. 

d Angle from the center of the SNR. 

G Hardncss ratio, computed according to (Ch — Cl) / {Ch + Cl), where C is the number of counts in a band. 



Note. — Quantities in parentheses are 1-cr uncertainties. 
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Table 6 

X-ray Sources in SNR G084.2-0.8 



ID" 


CXOU J b 




a 


(J2000) 


6 






rgo" 
(") 


AR d 

(') 


Counts L 


Counts H 


CountSA 


Counts 
(0.5 - 2.0 koV) 


HRl.h" 


1 


205328.9+432659 


20 


:53:28 


.96 


+43: 


26 


:59 


.1 


0.5 


1.3 


40(6) 


-0.3(0.1) 


40(6) 


38(6) 


-1.01(0.01) 


2 


205357.9+432459 


20 


:53:57 


.92 


+43 


■24 


:59 


3 


1.3 


6.9 


5.9(2.8) 


37(6) 


44(7) 


7.2(3) 


0.73(0.12) 


3 


205349.1+432551 


20 


:53:49 


.18 


+43 


:25 


:51. 


.2 


0.9 


5.1 


11(3) 


22(5) 


33(6) 


11(3) 


0.33(0.17) 


4 


205250.5+433030 


20 


:52:50 


.55 


+43 


■30 


:30. 


.0 


1.2 


6.5 


27(5] 


4 5f2 4) 


31(6) 


27(5) 


-0.71(0.14) 


5 


205352.0+422809 


20 


:53:52 


.03 


+43 


:28 


:09 


.7 


1.2 


5.5 


6.1(2.6) 


14(4) 


20(5) 


5.8(2.6) 


0.40(0.21) 


6 


205306.9+423259 


20 


:53:06 


.89 


+43 


32 


:59 


.7 


1.2 


6.3 


-0.5(0.2) 


27(6) 


26(6) 


-0.5(0.2) 


1.04(0.02) 


7 


205248.2+433214 


20 


:52:48 


.25 


+43' 


.32 


:14 


.0 


2.0 


7.9 


4.3(2.4) 


23(5) 


26(6) 


4.0(2.4) 


0.69(0.16) 


10 


205351.7+432537 


20 


:53:51 


.76 


+43 


.25 


:37. 


.9 


1.4 


5.7 


1.7(1.8) 


14(4) 


16(4) 


1.7(1.8) 


0.79(0.21) 


11 


205317.6+432206 


20 


:53:17 


.66 


+43 


:22 


:06 


.8 


1.2 


5.2 


0.3(1.1) 


13(4) 


13(4) 


0.3(1.1) 


0.95(0.16) 


12 


205353.1+432751 


20 


:53:53 


.13 


+43 


.27 


:51 


.5 


1.1 


5.7 


8.7(3.2) 


4.0(2.4) 


13(4) 


8.7(3.2) 


-0.37(0.30) 


14 


205302.8+433207 


20 


:53:02 


.87 


+43 


.32 


:07 


.6 


1.5 


5.9 


12(4) 


2.4(1.8) 


13(4) 


12(4) 


-0.67(0.23) 


15 


205306.5+432822 


20 


53:06 


.50 


+43 


28 


:22 


.2 


0.6 


3.0 


4.2(2.1) 


11(3) 


16(4) 


4.2(2.1) 


0.46(0.23) 


16 


205335.6+433427 


20 


:53:35 


.66 


+43' 


:34 


:27 


.6 


1.9 


7.6 


14(4) 


0.4(1.9) 


12(4) 


9.3(3.5) 


-0.94(0.26) 


18 


205346.0+433345 


20 


:53:46 


.04 


+43: 


33 


:45 


.7 


2.4 


7.8 


-0.4(1.5) 


20(5) 


20(5) 


-0.1(1.5) 


1.04(0.16) 


22 


205242.6+432451 


20 


:52:42 


61 


+43 


:24 


:51 


.0 


1.2 


7.6 


9.9(3.7) 


4.9(2.8) 


15(5) 


8.3(3.3) 


-0.34(0.30) 


23 


205332.2+432355 


20 


:53:32 


.27 


+43' 


:23 


:55 


.7 


0.7 


3.9 


-0.3(0.1) 


17(4) 


17(4) 


-0.3(0.1) 


1.03(0.02) 



a Intcrnal identifier of the form SNR. G084.2-0.8: N. 

b Official IAU name. 

c Approximatc 90% confidence radius. 

Angle from the center of the SNR. 
°Hardncss ratio, computed according to (Ch — Cl)/(Ch + Cl), where C is the number of counts in a band. 



Note. — Quantities in parentheses arc 1-cr uncertainties. 



Table 7 

X-ray Sources in SNR G127.1+0.5 



ID a 


CXOU J b 


a 8 
(.72000) 


rga" 
(") 


AR d 

(') 


CountSL 


CountSH 


CountSA 


Counts 0. 5-2.0 koV 


HRl,h° 


1 


012830.6+630629 


01:28:30.64 


+63:06:29.9 


0.4 


0.2 


205(15) 


580(20) 


780(30) 


182(14) 


0.48(0.03) 


3 


012736.6+630345 


01:27:36.60 


+63:03:45.5 


1.3 


6.9 


27(5) 


40(7) 


68(9) 


25(5) 


0.19(0.13) 


4 


012807.7+630150 


01:28:07.72 


+63:01:50.5 


1.3 


5.5 


16(4) 


25(5) 


41(7) 


17(4) 


0.21(0.16) 


5 


012842.3+630825 


01:28:42.31 


+63:08:25.6 


0.6 


2.2 


9.5(3.2) 


21(5) 


32(6) 


8.4(3) 


0.38(0.17) 


8 


012925.4+630535 


01:29:25.43 


+63:05:35.3 


1.7 


6.1 


11(4) 


15(4) 


26(5) 


9.3(3.3) 


0.17(0.21) 


12 


012910.8+631014 


01:29:10.86 


+63:10:14.8 


1.6 


5.7 


10(3) 


8.5(3.2) 


20(5) 


11(3) 


-0.10(0.25) 


13 


012813.9+630621 


01:28:13.91 


+63:06:21.1 


0.9 


2.1 


5.3(2.4) 


11(3) 


17(4) 


6.3(2.6) 


0.33(0.24) 


14 


012755.2+630141 


01:27:55.23 


+63:01:41.6 


1.5 


6.4 


12(4) 


19 5 


33(6) 


7.9(3) 


0.23(0.18) 


15 


012837.8+630602 


01:28:37.83 


+63:06:02.8 


0.8 


0.8 


2.9(1.8) 


10(3) 


13(4) 


2.9(1.8) 


0.56(0.24) 


17 


012735.7+630241 


01:27:35.75 


+63:02:41.6 


1.7 


7.4 


2.0(2.1) 


17(4) 


19(5) 


0.1(1.5) 


0.79(0.21) 


18 


012917.4+630242 


01:29:17.42 


+63:02:42.3 


2.1 


6.4 


9.5(3.3) 


10(4) 


20(5) 


7.7(3) 


0.03(0.25) 


22 


012801.8+631244 


01:28:01.83 


+63:12:44.5 


2.0 


7.1 


6.4(2.8) 


5.9(3) 


12(4) 


5.3(2.6) 


-0.04(0.33) 



internal identifier of the form SNR G127.1+0.5:iV. 

b Official IAU name. 

c Approximatc 90% confidence radius. 

d Angle from the center of the SNR. 

°Hardncss ratio, computed according to (Ch — Cl) / (Ch + Cl), where C is the number of counts in a band. 



Note. — Quantities in parentheses arc 1-a uncertainties. 
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majority are not. These are meant to be internal designa- 
tions only and do not replace the official IAU names of the 
form CXOU JHHMMSS.siDDMMSS, where HHMMSS.s 
represents the Right Ascension and ±DDMMSS represents 
the Declination. When detailing the source identifications 
(Tabs. 4-7) we give both the internal and official names, 
but in the rest of this paper we use only the internal names. 
When initially identifying X-ray sources we numbered then 
consecutively, but the numbers in Tables 4-7 are no longer 
consecutive as we have removed sources at radii > 8' and 
with fewer than 10 counts (§ 4.2). 

4.3. X-ray Extended- Source Detection 

No extended sources were detected during the wavdetect 
runs, with a maximum wavelet scale of 16 pixels or 8". We 
also manually examined the X-ray images for sources with 
larger sizes and did not find any (there is some extended 
emission in SNR G315.4— 2.3 toward the south and west, 
but that is almost certainly due to diffuse thermal emis- 
sion from the SNR/RCW 86 complex as it has the same 
general spectrum; see Fig. 3). 

To quantify the limits on extended emission, we deter- 
mine the average background counts in a region free from 
sources. These counts are presented in Table 8, where we 
also find the 3 a limits to extended emission. 

4.4. X-ray Timing 

The majority of the detected sources have too few pho- 
tons for meaningful analyses of their variability. We can 
perform some analysis, though, with the brightest sources 
(we set a limit of 100 counts for a lightcurve, which 




,3 



10° 10 1 10 2 10 3 



Counts (0.3-2.0 keV) 

Fig. 10.— L-band count vs. H-band counts for SNR G093. 3+6.9. 
The sources from Table 4 are shown in green and are numbered. 
Model spectra (computed using PIMMS) with Nn = 2 X 10 21 cm~ 2 
are plotted as diagonal lines: power-law models as red dash-dotted 
lines (r is as indicated) and blackbody models as blue dashed lines 
(fcToo is as indicated). The solid thick black line represents the 
median spectrum from the CDF/Orion studies, with the dashed 
thick black lines showing the 25- and 75-percentile spectra. 




Counts (0.3-2.0 keV) 



Fig. 11 — L-band count vs. H-band counts for SNR G315.4-2.3. 
The sources from Table 5 are shown in green and are numbered. 
Model spectra (computed using PIMMS) with JVjj = 2 X 10 21 cm -2 
are plotted as diagonal lines: power-law models as red dash-dotted 
lines (r is as indicated) and blackbody models as blue dashed lines 
(fcToo is as indicated). The solid thick black line represents the 
median spectrum from the CDF/Orion studies, with the dashed 
thick black lines showing the 25- and 75-percentile spectra. 




Counts (0.3-2.0 keV) 



Fig. 12.— L-band count vs. H-band counts for SNR G084.2-0.8. 
The sources from Table 6 are shown in green and are numbered. 
Model spectra (computed using PIMMS) with JVh = 2 x 10 21 cm" 2 
are plotted as diagonal lines: power-law models as red dash-dotted 
lines (r is as indicated) and blackbody models as blue dashed lines 
(fcToo is as indicated). The solid thick black line represents the 
median spectrum from the CDF/Orion studies, with the dashed 
thick black lines showing the 25- and 75-percentile spectra. 



A Search for Compact Central Sources in Four SNRs 15 



Table 8 

Limits to Extended X-ray Emission in SNRs G093.3+6.9, G315.4-2.3, G084.2+0.8, & G127. 1+0.5 



SNR 


BG Counts 


Count Limits 


L Limits 


L(6 = V) Limits 




(xlO~ 2 arcsec -2 ) 


(xlO- 2 ) 


(x6» 2 10 28 ergs s" 1 ) 


(xlO 32 ergs s~ x ) 


G093.3+6.9 


1.75(1) 


66> 2 + 700 


8 


3 


G315.4-2.3 


2.10(2) 


76> 2 + 700 


10 


4 


G084.2+0.8 


2.34(2) 


66> 2 + 706» 


10 


4 


G127.1+0.5 


2.03(2) 


76» 2 + 806» 


1 


0.3 



Note. — 9 is the radius of the extended region in arcseconds. The limits assume uniform 
weighting over 9, and are at the 3cr level. The counts are in the 0.3-8.0 keV range. The 
luminosity limits assume a power-law with T — 1.5. 




Counts (0.3-2.0 keV) 



Fig. 13.— L-band count vs. H-band counts for SNR G127. 1+0.5. 
The sources from Table 7 are shown in green and are numbered. 
Model spectra (computed using PIMMS) with ./Vh = 2 x 10 21 cm~ 2 
are plotted as diagonal lines: power-law models as red dash-dotted 
lines (r is as indicated) and blackbody models as blue dashed lines 
(fcToo is as indicated). The solid thick black line represents the 
median spectrum from the CDF/Orion studies, with the dashed 
thick black lines showing the 25- and 75-percentile spectra. 



eliminated all of the sources in SNRs G315.4— 2.3 and 
G084.2+0.8). In Figure 14 we show X-ray lightcurves for 
the two sources in SNR G093.3+6.9 with > 100 counts, 
SNR G093.3+6.9:l and SNR G093.3+6.9:2. The lightcurve 
for SNR G093.3+6.9:2 is quite constant, but that for 
SNR G093.3+6.9:l has a significant flare lasting w 1.5 hr. 
Searches for periodic variation, however, showed nothing. 
With the 3.2-s sampling of ACIS-I we were unable to 
search for rapid variability, but the low count-rates made 
that impossible anyway. We took the barycentered arrival- 
time data for SNR G093.3+6.9:l and SNR G093.3+6.9:2 
and performed FFTs with the data binned into 4-s inter- 
vals. No significant peaks in the periodogram were found 
over the frequency range 0.01-0.12 Hz. 

In Figure 15 we show the X-ray lightcurve for the 
only source in SNR G127. 1+0.5 that has > 100 counts: 
SNR G127.1+0.5:l. The lightcurve is consistent with a 
constant flux. 

4.5. Multi-wavelength Followup & Counterpart 
Identification 

4.5.1. Motivation 

After identifying X-ray sources with Chandra (or 
XMM), the question is then to determine which, if any, 
are the compact remnants of the SNRs. We have used 
source-count statistics (Brandt et al. 2001; Mushotzky 
et al. 2000; Feigelson et al. 2002) to estimate the number 




Fig. 14. — Background-subtracted X-ray lightcurves of sources 
in SNR G093.3+6.9 with > 100 counts. Left: SNR G093.3+6.9:l; 
right: SNR G093.3+6.9:2. The dashed red line is 0.3-2.0 keV (L- 
band), the dot-dashed blue line is 2.0-8.0 keV (H-band) and the 
solid black line is 0.3-8.0 keV (A-band). A flare lasting » 6000 s 
and with an amplitude change of ss 500% is present in the lightcurve 
of SNR G093.3+6.9:!. 
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of foreground/background sources given the A^h and dif- 
fuse SNR background toward each target — we expect to 
detect 10-50 field sources toward each SNR in the CXO- 
A03 sample. Most of these sources will be detected with 
few counts (10-30), and will therefore not be suitable for 
X-ray spectroscopy that could, of itself, determine which 
are neutron stars. Because of their small X-ray count- 
rates, weeding out interlopers requires multi-wavelength 
observations. 

Isolated neutron stars have high X-ray to optical flux ra- 
tios (Hulleman, van Kerkwijk, & Kulkarni 2000) — deep 
optical/IR imaging is therefore an efficient way to identify 
background sources. We thus follow our X-ray observa- 
tions with successively deeper optical and IR observations, 
identifying progressively fainter counterparts as we go. 

Interloper sources, on the other hand, typically have 
much brighter optical/IR counterparts (Hertz & Grindlay 
1988; Grindlay et al. 2003). In the Galactic plane, the 
majority of sources are either nearby bright stars or active 
late-type stars. Other possible sources are RS CVn bina- 
ries, X-ray binaries, or cataclysmic variables. The extra- 
galactic sources are usually AGN or star-forming galaxies, 
although some nearby spiral galaxies are also detectable 
(e.g., Barger et al. 2003). Stars with spectral types A- 
M can have some detectable X-ray emission (largely de- 
pendent on rotation and magnetic fields) that appears as 
a hot (5-7 x 10 6 K) plasma (Katsova & Cherepashchuk 
2000). X-ray binaries usually have harder spectra char- 
acterized by power-laws with indices T rj 1.5-2 (Grindlay 
et al. 2003). Stars are unresolved optical/IR sources, but 
for nearby binaries the X-ray emission may not lie exactly 
coincident with the optical/IR emission (if, for example, 
the optical emission is from the stars but the X-ray emis- 
sion is from interacting coronae between the stars) and the 
binary members may be merged to give the appearance of 
an extended sources. The extragalactic sources also have 
hard spectra with power-law indices T « 1.5-2 (Bauer 
et al. 2002), and are optically fainter (R > 18). While the 
nuclei of these galaxies are unresolved, the closer galaxies 



0.06 



0.05 



0.04 



0.01 




Fig. 15. — Background-subtracted X-ray lightcurve of 
SNR G127.1+0.5:l. The dashed red line is 0.3-2.0 keV (L-band), 
the dot-dashed blue line is 2.0—8.0 keV (H-band) and the solid black 
line is 0.3-8.0 keV (A-band). 



can have resolved optical/IR emission. 

There are a small number of galaxies that have some- 
what extreme X-ray-to-optical flux ratios. In a sample 
of 503 X-ray sources over 448 arcmin 2 , Koekemoer et al. 
(2004) have identified a few sources with X-ray-to-optical 
(850 nm) flux ratios exceeding 100, including 7 sources 
without 850 nm detections. While the identities are un- 
clear, they suggest that these galaxies have such faint op- 
tical counterparts through a combination of intrinsic red- 
dening and high redshift. However, these galaxies should 
not significantly impact our sample. We might expect 1-3 
of these sources in each of our Chandra images, but even 
without optical detections we should be able to identify 
these as galaxies in the near-infrared, as all but one of the 
sources were detected at K s band (where we have our deep- 
est observations) and the X-ray spectra are harder than 
those of most neutron stars (also see Gandhi et al. 2003). 
In addition, these galaxies appear to just be the tail end 
of the X-ray-to-optical flux ratio distribution, with most 
neutron stars being significantly higher. 

4.5.2. Execution 

Therefore, after the detection of the X-ray sources we 
obtained progressively deeper optical and IR data of the 
SNR fields. Starting with small telescopes (Palomar 60- 
inch, Las Campanas 40- inch) we moved to larger telescopes 
(Palomar 200-inch, NTT, Magellan, and finally Keck) as 
we identified counterparts to more and more of the X-ray 
sources (see below). The details of the optical/IR obser- 
vations arc in Tables 9-12. 

We began by registering all of the data for a SNR (we 
also included 2MASS J- and fC s -band images) to the same 
frame, and then searched for optical/IR counterparts to 
the X-ray sources using sextractor (Bertin & Arnouts 
1996). We considered a source as a potential counterpart 
if it was within the 90% confidence radius of the X-ray 
source combined with a 0'.'2 uncertainty for the optical as- 
trometry. We also inspected all X-ray sources to see if 
there was an optical source at the edge of the error cir- 
cle, if there were multiple sources, etc. We then deter- 
mined the positions and magnitudes of all of the detected 
sources, which we list in Table 13-16. The positions in 
Tables 13-16 are the averages of the positions for all im- 
ages where the source was detected, except for sources 
that were saturated in several bands but which had un- 
ambiguous 2MASS detections, where we used the 2MASS 
position alone. For sources that had multiple detections in 
the same band but by different instruments (e.g. C40 and 
EMMI i?-band) , we used the detection that had the high- 
est S/N without being saturated. We show postage-stamp 
images of the optical/IR counterparts to the X-ray sources 
in Figures 16-18 (for SNR G093. 3+6.9), Figures 19-21 (for 
SNR G315.4-2.3), Figures 22-24 (for SNR G084.2-0.8), 
and Figures 25-26 (for SNR G127. 1+0.5). The instru- 
ment(s)/band(s) chosen for each source are those that best 
illustrate the detection. 

To determine limiting magnitudes, we simulated Gaus- 
sian stars with the same FWHM as the average seeing for 
the image. We placed the stars randomly in regions that 
were not too crowded, similar to the regions where the X- 
ray sources actually were, and determined the 3 a limiting 
flux by performing photometry (again using sextractor) 
on the stars. 
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Table 9 

Optical/IR Observations of SNR G093.3+6.9 



Telescope 


Instrument 


Band(s) 


Date 


Seeing 


Exposure 


Sources' 1 








(UT) 


(arcsec) 


(sec) 




P60 


P60CCD 


BVRI 


2001-Jul-23 


1.9-2.5 


300 


-4 


P200 


LFC 


g'r' 


2002-Jan-18 


1.9 


600 


-13 


Keck I 


LRIS 


9 


2002-Jun-15 


0.9 


2300 


5,14,15,17 






R 




0.9 


2000 


5,14,15,17 






9 




0.9 


1380 


2,5,7,8,9,16,19,26 






R 




1.0 


1200 


2,5,7,9,16,19,26 


Keck I 


NIRC 


K s 


2002-Jun-02 


0.5 


250 


5,6,9,14,15,17 


P60 


P60IR 


K s 


2002-Jul-26 


2.0 


900 


4,10,13,19 


Keck II 


ESI 


R 


2002-Sep-03 


0.8 


540 


4 






I 




0.8 


900 


4 






I 




0.8 


2520 


8,26 






R 




0.7 


660 


10,24 






I 




0.6 


600 


10 






R 




0.7 


960 


13 


Keck I 


NIRC 


K s 


2002-Nov-16 


0.5 


2950 


19 


P200 


WIRC 


K s 


2003-Jul-24 


0.8 


1080 


-6,-10,-16 



a Indicatcs which sources from Table 4 were on which images. Negative numbers indicate that all 
sources but the negated one(s) were on the image. 

Note. — The telescopes/instruments used were P60CCD: the CCD camera on the Palomar 60- inch; 
P60IR: the Infrared camera on the Palomar 60-inch (Murphy et al. 1995); LFC: the Large Format 
Camera on the Palomar 200- inch; LRIS: the Low- Resolution Imaging Spectrometer on the 10-m 
Keck I telescope (Oke et al. 1995); NIRC: the Near-Infrared Camera on the 10-m Keck I telescope 
(Matthews & Soifer 1994); ESI: the Echellette Spectrograph and Imager on thelO-m Keck II telescope 
(Sheinis et al. 2002); WIRC: the Wide-field Infrared Camera on the Palomar 200-inch. 
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Table 10 

Optical/IR Observations of SNR G315.4-2.3 



Telescope Instrument Band(s) 


Date 


Seeing Exposure Sources 3, 




(UT) 


(arcsec) (sec) 



C40 


C40CCD 


BVRI 


2002-Apr-18 


1.3 


3600 


-3 






BR 


2002-Apr-20 


1.5 


3600 


-21 


NTT 


EMMI 


VRI 


2002-Jun-14 


0.8 


1130 


1,5,17,23 


Magellan II 


MagIC 


BR 


2003-Apr-03 


0.8 


1800 


1 






R 


2003-Apr-03 


0.8 


3000 


12,28 






R 


2003-Apr-03 


0.8 


3000 


6,17 






I 


2003-Apr-04 


0.7 


1500 


1 






I 


2003-Apr-04 


0.6 


900 


6,17 


Magellan II 


PANIC 


K s 


2003-Apr-18 


0.8 


2160 


1 






K s 


2003-Apr-18 


0.9 


2160 


12,28 






K s 


2003-Apr-18 


1.0 


2160 


6,17 






K s 


2003-Apr-20 


0.6 


2160 


1 






K s 


2003-Apr-20 


0.5 


2160 


12,28 






K a 


2003-Jun-07 


0.4 


900 


5 b 



Tndicates which sources from Table 5 were on which images. Negative numbers indicate 
that all sources but the negated one(s) were on the image. 

b Observed by M. van Kerkwijk. 



Note. — The telescopes/instruments used were C40CCD: the direct-imaging CCD camera 
on the Las Campanas 40-inch; EMMI: the ESO Multi-Mode Instrument (red imaging arm 
only) on the 3.5-m New Technology Telescope (NTT) at La Silla; MagIC: Raymond and 
Beverly Sackler Magellan Instant Camera on the 6.5-m Clay (Magellan II) telescope; PANIC: 
Persson's Auxiliary Nasmyth Infrared Camera on the 6.5-m Clay (Magellan II) telescope. 



Table 11 

Optical/IR Observations of SNR G084.2-0.8 



Telescope 


Instrument 


Band(s) 


Date 


Seeing 


Exposure 


Sources 3, 








(UT) 


(arcsec) 


(sec) 




P200 


LFC 


r' 


2002-Oct-07 


1.5 


2730 


all 


P200 


WIRC 


J 


2002-Nov-Ol 


1.2 


4500 


all 






K s 


2002-Nov-Ol 


0.9 


2640 


all 


Keck I 


NIRC 


K s 


2003- Aug- 11 


0.5 


450 


2,6,10,11,14,18,23 



indicates which sources from Table 6 were on which images. Negative numbers indicate that all 
sources but the negated one(s) were on the image. 



Note. — The telescopes/instruments used follow from Table 9. 
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Table 12 

Optical/IR Observations of SNR G127. 1+0.5 



Telescope 


Instrument 


Band(s) 


Date 


Seeing 


Exposure 


Sources 3, 








(UT) 


(arcsec) 


(sec) 




P60 


P60CCD 


R 


2001-Jul-23 


1.9 


3900 


-22 






1 


2001-Jul-23 


1.5 


300 


-22 


P200 


LFC 


r 


2002-Oct-07 


1.9 


2700 


all 


P200 


WIRC 


J 


2002-Oct-28 


0.9 


4050 


all 






K s 


2002-Oct-28 


0.8 


2700 


all 


Keck I 


NIRC 


K s 


2003- Aug- 11 


0.6 


900 


3,5,13,22 



indicates which sources from Table 7 were on which images. Negative numbers indicate 
that all sources but the negated one(s) were on the image. 



Note. — The telescopes/instruments used follow from Table 9. 



We define optical/IR fluxes using the zero-point cal- 
ibrations of Bessell, Castelli, & Plez (1998), where 
F(m) = F l(T m / 2 - 5 andF = (2.0 x 1(T 5 , 1.4 x 1CT 5 , 9.0 x 
10~ 6 ,8.7 x 10~ 7 ) ergs s" 1 cm" 2 for (V,R,I,K S ), respec- 
tively. We also use the reddening coefficients of Bessell 
et al. (1998), such that A x /A v = (1.0,0.82,0.62,0.11) 
again for (V, R, I, K s ), respectively. 

4.5.3. Counterpart Evaluation 

It is certainly possible that there will be an unrelated op- 
tical/IR source in the error circle of an X-ray source (e.g., 
Kaplan, Kulkarni, & Murray 2001), especially if the X-ray 
sources are off-axis and/or have few counts, and therefore 
have large position uncertainties. There is no entirely ac- 
curate way to prevent this from happening and leading 
to the false association of what is actually a neutron star 
with another source, which would cause us to reject the 
neutron star from our sample. 

However, there are ways that we can guard against 
this and incorporate our uncertainties into our limits. If 
enough is known about the optical/IR source (e.g. col- 
ors), it is possible to determine what the source is (Ka- 
plan et al. 2001) and thereby assign some likelihood to its 
association with an X-ray source. This is not always pos- 
sible for the sources in our sample as the information is 
often limited to optical/IR detections in one or two bands 
(Tabs. 13- 16). Similarly, for the few sources with enough 
X-ray counts a determination can be made based on the X- 
ray spectrum (or just hardness) and lightcurve (Fig. 14). 
For the majority of the sources, though, we must examine 
them relative to the sources in other samples, specifically 
the CDF and Orion samples (many of the stars in the 
Orion sample are younger than the general Galactic pop- 
ulation, but they do cover a wide range of stellar masses). 
For this reason we plot the optical/IR-to-X-ray flux ra- 
tio against X-ray flux in Figures 27-30 (we compute the 
X-ray flux from the observed count rate by using conver- 
sion factors determined from W3PIMMS for a blackbody 
with fcToo = 0.25 keV and Ah appropriate for each SNR; a 
change to a power-law with photon index of 1.5 raises the 
implied fluxes by a factor of m 2). Sources that fall in the 
loci defined by the other samples are likely to be of similar 



type. This is of course not a definitive assignment, but it 
should work most of the time. For instance, if we had as- 
sociated Star A from Kaplan et al. (2001) with the X-ray 
flux from CXO J232327.8+584842 (the central source in 
Cas A), it would appear in Figure 27 with the same X-ray 
flux as the compact source in Cas A and < 1 magnitude 
brighter than the limit we have plotted. In other words, it 
would be far outside the Orion/CDF loci. This argument, 
that Star A could not be the source of the X-rays from 
CXO J232327.8+584842 just on the basis of its X-ray-to- 
optical flux ratio, was made by Kaplan et al. (2001), who 
then followed it up by a detailed analysis of multi-band 
photometry of Star A (see also Fesen et al. 2002). 
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Fig. 16. — Optical/IR images of counterparts to X-ray sources in SNR G093.3+6.9. North is up, East is to the left, and a scale-bar is in 
the upper-left corner. 
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Fig. 17. — Optical/IR images of counterparts to X-ray sources in SNR G093.3+6.9 (cont.). North is up, East is to the left, and a scale-bar 
is in the upper-left corner. 
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Fig. 18. — Optical/IR images of counterparts to X-ray sources in SNR G093.3+6.9 (cont.). North is up, East is to the left, and a scale-bar 
is in the upper-left corner. 
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Fig. 19. — Optical/IR images of counterparts to X-ray sources in SNR G315.4— 2.3. North is up, East is to the left, and a scale-bar is in 
the upper-left corner. 
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Fig. 20. — Optical/IR images of counterparts to X-ray sources in SNR G315.4— 2.3 (cont.). North is up, East is to the left, and a scale-bar 
is in the upper-left corner. 
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Fig. 22. — Optical/IR images of counterparts to X-ray sources in SNR G084.2— 0.8. North is up, East is to the left, and a scale-bar is in 
the upper-left corner. 
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Fig. 23. — Optical/IR images of counterparts to X-ray sources in SNR G084.2— 0.8 (cont.). North is up, East is to the left, and a scale-bar 
is in the upper-left corner. 
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Fig. 24. — Optical/IR images of counterparts to X-ray sources in SNR G084.2— 0.8 (cont.). North is up, East is to the left, and a scale-bar 
is in the upper-left corner. 
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Fig. 25. — Optical/IR images of counterparts to X-ray sources in SNR G127. 1+0.5. North is up, East is to the left, and a scale-bar is in 
the upper-left corner. 
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Fig. 26. — Optical/IR images of counterparts to X-ray sources in SNR G127. 1+0.5 (cont.). North is up, East is to the left, and a scale-bar 
is in the upper-left corner. 
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Table 13 

Optical/IR Matches to X-ray Sources in SNR G093.3+6.9 



ID° 












5 b 






Ar c 


B 


v 


R 




g/g' d 




j 


K 


N(< K )° 










(.T2000) 








(arcscc) 








(mag) 








1 


20 


52 


: 22 


00 


+ 55: 


23: 


43. 


7 


0.1 




TO 1 c/o\ 

1 o. ±o{o) 


1 1 7c 1 n\ 
IZ. 1 o(o) 


12.535(10) 


sat 


sat 


ii.i ( 2^yj 


lU.ooy (1 1 ) 


< 0.01 


2 


20 


52 


:31. 


01 


+ 55: 


14 


:37. 


6 


0.3 


>23.5 


>22.5 


20.74(5) 


19.52(4) 


22.61(10) 


21.05(6) 


>15.8 


15.52(1) 


< 0.01 


4 


20 


52 


:25 


.83 


+ 55: 


27 


:41 


7 


0.1 






23.30(11) 


21.60(10) 


>24.0 


>23.5 


>15.8 


17.46(17) 


< 0.01 


5 


20: 


52 


:22 


01 


+ 55: 


:15 


:16. 


9 


0.1 


>23.5 


>22.5 


>25.0 


>20.7 


>26.5 


>23.5 


>15.8 


20.1(4) 


< 0.01 


6 


20: 


51 


:31. 


33 


+ 55: 


.20 


:31. 


4 


0.4 


>23.5 


>22.5 


>22.3 


>20.7 


>24.0 


>23.5 


>15.8 


16.1(4) 


< 0.01 


7 


20 


52 


:50. 


22 


+ 55: 


16 


:06. 


7 


0.2 


17.67(7) 


16.23(3) 


15.18(3) 


13.857(10) 


sat 


sat 


12.363(16) 


11.475(16) 


< 0.01 


8 


20 


52 


:42 


13 


+ 55: 


14 


:10. 


7 


1.1 


>23.5 


>22.5 


>22.3 


23.43(13) 


>26.1 


>23.5 


>15.8 


18.9(1) 


0.06 


9 


20 


52 


:36 


.48 


+55. 


:17 


:45. 


8 


0.2 


>23.5 


>22.5 


>25.1 


>20.7 


>26.1 


>25.1 


>15.8 


22.3(4) 


< 0.01 


10 


20 


51 


:39 


.85 


+ 55: 


■25 


:53. 


1 


0.3 


>23.5 


>22.5 


22.99(11) 


21.36(10) 


>24.0 


>23.5 


>15.8 


17.45(11) 


< 0.01 


13 


20: 


52 


:42 


.23 


+55. 


26 


:07. 


4 


0.4 






22.5(5) 








>15.8 


17.44(6) 


< 0.01 


14 


20 


52 


:05 


08 


+ 55: 


:17 


:58. 


6 


0.2 


>23.5 


>22.5 


22.72(6) 


>20.7 


>26.3 


>23.5 


>15.8 


17.04(5) 


< 0.01 


15 


20 


52 


:09 


.83 


+ 55: 


:18 


:21 


4 


0.4 


>23.5 


>22.5 


>22.3 


>20.7 


>26.3 


>25.2 


>15.8 


20.2(4) 


0.01 


16 


20 


52 


:50. 


.16 


+ 55: 


20 


:25. 


3 


0.7 


17.63(7) 


16.22(3) 


15.49(3) 


14.732(10) 


sat 


sat 


13.61(3) 


13.02(3) 


< 0.01 


17 


20 


52 


:12. 


66 


+55. 


18 


:54. 


6 


0.2 


>23.5 


>22.5 


>25.2 


>20.7 


>26.3 


>23.5 


>15.8 


21.1(4) 


< 0.01 


19 


20 


52 


:36 


.79 


+ 55: 


15 


:28. 


1 


1.1 


>23.5 


>22.5 


24.99(14) 


>20.7 


>26.1 


>23.5 


>15.8 


19.9(4) 


0.09 




20: 


52 


:36 


02 


+55. 


:15 


:30 


.5 


1.8 


>23.5 


>22.5 


>22.3 


>20.7 


>26.1 


>23.5 


>15.8 


19.5(4) 


0.22 


24 


20 


51 


:52. 


31 


+ 55: 


26: 


:02. 


8 


0.7 


>23.5 


>22.5 


25.2(2) 


>20.7 


>24.0 


>23.5 


>15.8 


19.1(1) 


0.03 


26 


20 


52 


:48 


.31 


+ 55: 


:14 


:21 


8 


0.5 


>23.5 


>22.5 


>22.3 


23.00(11) 


>26.1 


>23.5 


>15.8 


17.98(7) 


0.01 


30 


20 


51 


:29 


.88 


+55. 


15 


:48. 


2 


1.4 


>23.5 


>22.5 


>22.3 


>20.7 


>24.0 


>23.5 


>15.8 


18.1(1) 


0.07 



a Source number from Tabic 4. 

b Position of optical/IR source, averaged over all the bands in which there were detections. 

c Position difference between the X-ray and optical/IR source. 

For brevity, we give both (/-band (from LRIS) and (/'-band (from LFC) data in the same column. The limits of 24.0-mag are (/-band, while the 
limits of 26. 1/26.3-mag are (/-band. 

G The chance of finding a star within Ar of the X-ray source given the K s magnitude, using the star-count model of Nakajima et al. (2000). 



Table 14 

Optical/IR Matches to X-ray Sources in SNR G315.4-2.3 



ID a 


o 


P 






5 b 






Ar c 


B 


V 


R 




.7 


K„ 


7V(< K s ) d 








(.T2000) 








(arcsec) 






(mag) 






1 


14:43 


19 


.33 


-62 


:28 


:04. 


1 


0.2 


> 22.3 


> 24.5 


> 23.9 


> 22.5 


> 15.8 


19.60(8) 


0.04 


2 


14:43: 


33 


.73 


-62 


:29 


:27 


7 


0.5 


14.52(2) 


13.23(4) 


12.41(8) 


11.5(1) 


10.62(2) 


9.84(2) 


< 0.01 


3 


14:41: 


51 


.73 


-62 


:28 


:33. 


4 


1.0 


> 22.8 


19.88(6) 


> 15.8 


> 14.6 




5 


14:42: 


19 


.56 


-62 


:28 


:34. 


6 


0.2 


> 22.3 


> 24.5 


> 23.9 


> 22.5 


> 15.8 


20.5(4) 


0.05 


6 


14:43: 


:20 


.72 


-62 


:33 


:08. 


4 


0.7 


14.49(2) 


13.73(4) 


13.27(8) 


12.63(10) 


12.23(2) 


11.77(3) 


< 0.01 


8 


14:43- 


.40 


.99 


-62 


:31 


:38. 


3 


0.4 


sat 


12.39(4) 


11.92(8) 


11.35(10) 


10.88(2) 


10.50(2) 


< 0.01 


9 


14:43: 


46 


.51 


-62 


24 


:12. 


.1 


1.4 


20.11(3) 


18.58(3) 


17.35(4) 


16.12(7) 


14.50(6) 


13.40(5) 


0.02 


12 


14:43: 


:15 


.26 


-62 


:21 


:28. 


.7 


0.7 


> 22.3 


> 21.8 


> 25.5 


> 19.7 


> 15.8 


19.44(10) 


0.38 


13 


14:43: 


50 


.38 


-62 


:30 


:40. 


.8 


0.9 


19.71(3) 


17.79(5) 


17.04(8) 


15.55(10) 


14.33(5) 


13.31(5) 


0.01 


17 


14:43: 


21 


.90 


-62 


:32 


:10 


.0 


1.5 


> 22.3 


> 21.8 


> 25.6 


20.73(5) 


> 15.8 


17.71(6) 


0.43 


19 


14:42: 


10 


.79 


-62 


:22 


:02. 


9 


1.1 


13.63(2) 


13.10(2) 


12.92(4) 


12.29(7) 


11.52(2) 


11.19(3) 


< 0.01 


21 


14:42: 


52 


.36 


-62 


:21 


:08. 


.1 


1.7 


> 22.3 


> 21.8 


20.76(2) 


> 19.7 


> 15.8 


> 14.6 






14:42: 


52 


.10 


-62 


:21 


:06 


.5 


0.7 


> 22.3 


> 21.8 


19.69(8) 


> 19.7 


> 15.8 


> 14.6 




23 


14:42: 


13 


.11 


-62 


:32 


:20 


6 


0.2 


17.83(2) 


16.18(4) 


15.00(8) 


13.76(10) 


12.46(3) 


11.42(3) 


< 0.01 


26 


14:42: 


06 


.91 


-62 


:22 


:04. 


.7 


1.8 


17.18(2) 


16.26(2) 


15.53(4) 


14.87(7) 


14.02(4) 


13.37(6) 


0.03 


28 


14:43: 


20 


.25 


-62 


:21 


:38. 


.7 


1.3 


> 22.3 


> 21.8 


23.14(6) 


> 19.7 


> 15.8 


18.44(6) 


0.62 




14:43: 


19 


.88 


-62 


:21 


:40. 


.2 


2.0 


> 22.3 


> 21.8 


22.23(6) 


> 19.7 


> 15.8 


18.38(6) 


1.41 




14:43: 


.20 


.07 


-62 


:21 


:39. 


.4 


0.7 


> 22.3 


> 21.8 


18.66(5) 


> 19.7 


> 15.8 


16.12(5) 


0.02 



a Source number from Table 5. 

b Position of optical/IR source, averaged over all the bands in which there were detections. 

c Position difference between the X-ray and optical/IR source. 

The chance of finding a star within Ar of the X-ray source given the K s magnitude, using the star-count model of 
Nakajima ct al. (2000). 
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Table 15 

Optical/IR Matches to X-ray Sources in SNR G084.2-0.8 



ID a 




<5 b 


Ar c 


r' 


J 


K 


N(< K s ) d 




(J2000) 


(arcsec) 




(mag) 




1 


20:53:28.97 


+43:26:58.7 


0.4 


sat 


10.36(2) 


10.11(2) 


< 0.01 


2 


20:53:57.97 


+43:24:58.6 


0.9 


>23.2 


>20.8 


19.68(6) 


0.41 


3 


20:53:49.20 


+43:25:51.1 


0.2 


>23.2 


>20.8 


18.32(7) 


< 0.01 


4 


20:52:50.51 


+43:30:29.3 


0.8 


15.08(7) 


12.99(2) 


12.25(3) 


< 0.01 


5 


20:53:52.06 


+43:28:10.8 


1.2 


> 23.2 


17.44(4) 


15.97(2) 


0.05 




20:53:51.87 


+43:28:09.9 


1.7 


18.83(7) 


15.45(2) 


14.39(1) 


0.03 




20:53:51.99 


+43:28:08.9 


0.8 


>23.2 


>20.8 


17.95(5) 


0.10 


6 


20:53:06.86 


+43:32:59.3 


0.5 


>23.2 


>20.8 


19.45(6) 


0.11 


7 


20:52:48.32 


+43:32:12.7 


1.5 


>23.2 


20.67(13) 


17.64(3) 


0.27 


10 


20:53:51.79 


+43:25:37.1 


0.9 


>23.2 


20.69(15) 


17.10(4) 


0.06 




20:53:51.65 


+43:25:37.2 


1.3 


>23.2 


>20.8 


21.0(2) 


1.53 


11 


20:53:17.76 


+43:22:06.6 


1.1 


>23.2 


18.43(4) 


15.332(10) 


0.03 


12 


20:53:53.17 


+43:27:52.0 


0.7 


sat 


13.39(3) 


11.88(2) 


< 0.01 


14 


20:53:02.75 


+43:32:08.5 


1.5 


sat 


11.97(2) 


11.61(2) 


< 0.01 


15 


20:53:06.49 


+43:28:21.9 


0.3 


>23.2 


>20.8 


18.37(5) 


0.02 


16 


20:53:35.61 


+43:34:27.5 


0.6 


18.22(7) 


14.84(4) 


14.01(7) 


< 0.01 


18 


20:53:46.08 


+43:33:46.4 


0.9 


>23.2 


>20.8 


20.91(14) 


0.71 




20:53:45.89 


+43:33:44.8 


1.8 


>23.2 


>20.8 


18.08(3) 


0.56 


22 


20:52:42.55 


+43:24:50.5 


0.8 


19.40(7) 


14.73(7) 


13.26(5) 


< 0.01 


23 


20:53:32.25 


+43:23:55.2 


0.5 


>23.2 


>20.8 


18.23(4) 


0.05 



a Source number from Tabic 6. 

b Position of optical/IR source, averaged over all the bands in which there were detec- 
tions. 

c Position difference between the X-ray and optical/IR source. 

d The chance of finding a star within Ar of the X-ray source given the K s magnitude, 
using the star-count model of Nakajima et al. (2000). 
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Table 16 

Optical/IR Matches to X-ray Sources in SNR G127.1+0.5 



ID a 










5 h 






Ar e 


R 




r 


J 




N(< K s ) d 








(.T2000) 








(arcscc) 






(mag) 






1 


01 


:28:30. 


.58 


+63: 


06 


:29 


,7 


0.4 


16.59(4) 


sat 


17.76(5) 


13.96(12) 


12.09(8) 


< 0.01 


3 


01 


:27:36 


51 


+ 63: 


:03 


:45 


.5 


0.6 


>23.0 


>21.3 


>24.5 


blended 


17.82(3) 


0.03 


4 


01 


:28:07 


.70 


+ 63: 


01 


:50 


.6 


0.2 


22.18(9) 


20.69(10) 


22.64(7) 


19.45(6) 


17.74(4) 


< 0.01 


5 


01 


:28:42. 


25 


+ 63: 


08 


:25 


.8 


0.4 


>23.0 


>21.3 


>24.5 


>21.1 


18.80(4) 


0.02 


8 


01 


:29:25 


.42 


+ 63: 


05 


:34 


.2 


1.1 


21.80(8) 


20.48(8) 


22.43(6) 


18.39(2) 


16.86(3) 


0.04 


12 


01 


:29:10. 


.80 


+ 63: 


10: 


14 


.4 


0.6 


>23.0 


>21.3 


>24.5 


18.91(5) 


16.60(3) 


0.01 


13 


01 


:28:13. 


.86 


+ 63: 


06 


21 


.2 


0.3 


>23.0 


>21.3 


>24.5 


>21.1 


19.85(8) 


0.02 


14 


01 


:27:55 


.35 


+ 63: 


01 


:41 


.5 


0.8 


>23.0 


>21.3 


>24.5 


20.32(7) 


17.08(3) 


0.03 


15 


01 


:28:37 


.77 


+63. 


06 


03 


.3 


0.7 


>23.0 


>21.3 


>24.5 


blended 


18.53(8) 


0.06 


17 


01 


:27:35. 


.74 


+ 63: 


02 


:42 


,4 


0.8 


>23.0 


>21.3 


>24.5 


>21.1 


18.17(8) 


0.06 


18 


01 


:29:17. 


.27 


+ 63: 


:02 


:42 


.3 


1.0 


>23.0 


>21.3 


>24.5 


16.717(14) 


14.774(15) 


0.01 


22 


01 


:28:01 


81 


+63. 


12 


:44 


2 


0.3 


>23.0 


>21.3 


>24.5 


>21.1 


20.34(13) 


0.03 




01 


:28:01. 


93 


+ 63: 


12 


:45 


.7 


1.4 


>23.0 


>21.3 


>24.5 


>21.1 


20.81(15) 


0.75 



a Source number from Table 7. 

b Position of optical/IR source, averaged over all the bands in which there were detections. 
c Position difference between the X-ray and optical/IR source. 

d The chance of finding a star within Ar of the X-ray source given the K s magnitude, using the star-count 
model of Nakajima et al. (2000). 
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This will not work so well for the sources here, 
as the X-ray fluxes are significantly less than that of 
CXO J232327.8+584842 and consequently the location 
in Figure 27 for a random optical/IR counterpart would 
be closer to the Orion/CDF loci. 

To determine the probabilities of random false associa- 
tions, we have used the IR star-counts as modeled by Naka- 
jima et al. (2000) and galaxy-counts observed by Cimatti 
et al. (2002). We plot in Figure 31 the number of sources 
brighter than a given K s magnitude per square arcsecond 
for all four SNRs (due to their lower Galactic latitudes, the 
numbers for SNR G315.4-2.3 and SNR G084.2-0.8 are 
a factor of 5-15 higher than those for SNR G093.3+6.9, 
and SNR G127. 1+0.5 is in between). We chose to plot the 
counts in the K a band as these are least affected by extinc- 
tion. In Tables 13-16 we also give, for sources with K s de- 
tections, the expected number of random stars 13 brighter 
than the detected source within Ar (the distance between 
the K s source and the X-ray source). For most cases, es- 
pecially the bright 2MASS stars, these numbers are very 
low (< 0.01), but for a few of the fainter sources the num- 
bers can become significant (> 0.5). For those cases with 
roughly even chances of random associations we also exam- 
ine (in §§ 5.1.3, 5.2.3, 5.3.3, and 5.4.3) other factors like the 
source spectrum and optical/IR colors, where available. In 
Figures 32-34 we show plots of expected stellar color vs. 
magnitude for different distances and extinctions. Using 
these plots can help us determine the approximate type of 
the stellar counterparts in Tables 13-16, as demonstrated 
in Sections 5.1.3, 5.2.3, 5.3.3, and 5.4.3. These plots are 
not exact, though: they assume Ay cx D and they are 
only for main-sequence stars. When possible, the J — K s 
color should be used instead of other combinations as it is 
least susceptible to reddening effects (E(J — K s ) = 0.18, 
compared to E(R — K s ) = 0.71), but it is less sensitive to 
stellar type than other combinations. 

In Figures 35 and 36 we examine the X-ray-to-optical 
offset. We see that even for sources with large offsets, the 
values are reasonably consistent with the expected distri- 
bution. Compare to Figure 1 of Barger et al. (2003). 

We can see some general trends from the discussions 
for the different SNRs (SS 5.1.3, 5.2.3, 5.3.3, and 5.4.3) 
and from Figures 27-30. Namely SNR G093.3+6.9 and 
SNR G127.1+0.5 have the highest fractions of proba- 
ble galaxies, with SNR G084.2— 0.8 intermediate, and 
SNR G315.4-2.3 the lowest (summarized in Tab. 17). 
This is largely due to the differing Galactic latitudes and 
longitudes of the four SNRs, and is also seen somewhat 
in Figure 31. SNR G315.4-2.3, at a relatively low |6| 
and well in quadrant IV, is along the line of sight of 
many stars and the short exposure time of the observation 
helps to keep the number of background galaxies down. 
SNR G084.2-0.8 is at a lower but it's position near 
I = 90° and at the edge of a spiral arm (Feldt & Green 
1993) lowers the number of stars, while the longer exposure 
identified more galaxies (many of whom are heavily red- 
dened) . The total TVh along this like of sight helps to delin- 
eate galaxies and stars even without optical counterparts, 
as can be seen from Figure 12. Finally, SNR G093.3+6.9 

13 As seen in Figure 31, the chance of random associations with a 
galaxy is quite small compared to the chance of association with a 
star, and is usually negligible. 
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Fig. 31. — Number density of IR sources. The number per square 
arcsecond brighter than a given K B magnitude is plotted against 
K a magnitude. The data for stars are from the model of Naka- 
jima et al. (2000), and are plotted for SNR G093.3+6.9 (solid line), 
SNR G315.4-2.3 (dashed line), SNR G084.2-0.8 (dash-dotted line), 
and SNR G127.1+0.5 (dotted line). For galaxies (circles) the counts 
are from the K-20 project (Cimatti et al. 2002). The observed 
data for SNR G093.3+6.9, SNR G084.2-0.8, and SNR G127.1+0.5 
are from the WIRC observations, and are likely not complete for 
K s > 19. The predicted star counts are within a factor of 3 of the 
observed counts. 




Fig. 32. — Stellar color vs. magnitude for main-sequence stars, 
from Cox (2000, pp. 151 & 388): we have assumed that Ay oc d, 
normalized to Ay = 1 mag at d = 3 kpc (roughly appropriate for 
SNR G093.3+6.9 and SNR G315.4-2.3). Left: R — K 3 color vs. 
K s magnitude. Right: J — K s color vs. K s magnitude. Shown are 
tracks for luminosity class V stellar types F0, F5, GO, G5, K0, K5, 
M0, M2, and M5 (as labeled). Points along the tracks are labeled 
by (log 10 (ii/pc), Ay), and the tracks progress from (2.0,0.03) to 
(4.0,3.3). The maximum values of Ay expected along the lines 
of sight to SNR G093.3+6.9 and SNR G315.4-2.3 are 2 mag and 
5 mag, respectively (Predehl Sz Schmitt 1995). The majority of the 
stars from Nakajima et al. (2000) are of spectral type M 0-M3. 




Fig. 33. — Stellar color vs. magnitude for main-sequence stars, 
from Cox (2000, pp. 151 & 388): we have assumed that Ay oc d, 
normalized to Ay = 1 mag at d = 4.5 kpc (roughly appropriate 
for SNR G084.2— 0.8). Otherwise the figures are the same as Fig- 
ure 32. The maximum value of Ay expected along the line of sight 
to SNR G084.2-0.8 is 5 mag (Predehl & Schmitt 1995). 
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Fig. 27. — X-ray-to-optical/IR flux ratio vs. X-ray flux for sources in SNR G093.3+6.9 (Tabs. 4 and 13) with sources from the CDF/Orion 
studies and selected neutron stars. Upper left: V-band; upper right: i?-band; lower left: /-band; lower right: iC s -band. Stars from CDF/Orion 
are blue asterisks, galaxies are green circles. Selected neutron stars are black diamonds/limits, and are labeled. The unidentified X-ray sources 
in SNR G093.3+6.9 are the red squares/limits and are also labeled (in the case of multiple possible counterparts, the source is plotted multiple 
times). The diagonal lines represent constant magnitude, and are labeled by that magnitude. For the sources in SNR G093.3+6.9 the counts 
were converted to a flux by F 5 _ 2 q k c y = counts 5 _ 2 j^y X 3.4 X 10 — 16 ergs s _1 cm -1 , appropriate for a blackbody with fcToo = 0.25 keV 
and Nh = 2 X 10 21 cm -2 . As seen in Table 18, the X-ray fluxes change can be a factor of fa 2 higher for a power-law spectrum. 
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Fig. 28. — X-ray-to-optical/IR flux ratio vs. X-ray flux for sources in SNR G315.4— 2.3 (Tabs. 5 and 14) with sources from the CDF/Orion 
studies and selected neutron stars, following Figure 27. Upper left: V-band; upper right: i?-band; lower left: /-band; lower right: ii" s -band. 



For the sources in SNR G315.4— 2.3 the counts were converted to a flux by fo.5— 2.0 keV 
appropriate for a blackbody with kT x = 0.25 keV and A'jj = 2 X 10 21 cm -2 . 



countSQ 5 _2 o keV x 5.1 X 10 ergs s cm 
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Fig. 29. — X-ray-to-optical/IR flux ratio vs. X-ray flux for sources in SNR G084.2— 0.8 (Tabs. 6 and 15) with sources from the CDF/Orion 
studies and selected neutron stars, following Figure 27. Left: i?-band; right: if s -band. For the sources in SNR G093.3+6.9 the counts were 
converted to a flux by ^o. 5-2.0 keV = counts 5 _ 2 g keV x 2-0 X 10 — 16 ergs s _1 cm -1 , appropriate for a blackbody with kT x = 0.25 keV and 
N H = 2 X 10 21 cm" 2 . 
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Fig. 30. — X-ray-to-optical/IR flux ratio vs. X-ray flux for sources in SNR G127.1+0.5 (Tabs. 7 and 16) with sources from the CDF/Orion 
studies and selected neutron stars, following Figure 27. Left: H-band; right: iC s -band. For the sources in SNR G127. 1+0.5 the counts were 
converted to a flux by F 5 _ 2 g keV = countsg 5 _ 2 g k e y X 2.6 X 10 — 16 ergs s _1 cm -1 , appropriate for a blackbody with fcToo = 0.25 keV and 
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Fig. 34. — Stellar color vs. magnitude for main-sequence stars, 
from Cox (2000, pp. 151 & 388): we have assumed that Ay oc d, 
normalized to Ay = 1 mag at d = 1.3 kpc (roughly appropriate 
for SNR G127.1+0.5). Otherwise the figures are the same as Fig- 
ure 32. The maximum value of Ay expected along the line of sight 
to SNR G127.1+0.5 is 5 mag (Predehl & Schmitt 1995). 



Fig. 35. — Distribution of offset between X-ray sources and their 
proposed optical/IR counterparts as a function of off-axis angle. 
Left: offset (in arcsec) vs. off-axis angle (in arcmin). Right: offset 
normalized to the 90% radius (from Tabs. 4-7) vs. off-axis angle 
(in arcmin). Sources from SNR G093.3+6.9 are circles, those from 
SNR G315.4-2.3 are x's, those from SNR G084.2-0.8 are squares, 
and those from SNR G127.1+0.5 are stars. 
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Fig. 36. — Histogram of the distribution of offset between X-ray 
sources and their proposed optical/IR counterparts normalized to 
the 90% radius (from Tabs. 4- 7). Sources from SNR G093.3+6.9 
are dark blue, those from SNR G315.4— 2.3 are in white, those 
from SNR G084.2— 0.8 are the green hatched region, and those from 
SNR G127. 1+0.5 are the maroon hatched region. The red dashed 
line is the expected distribution (/(r) oc rexp(— r 2 )). 



is at a relatively high \b\ and also does not look toward the 
inner Galaxy, and SNR G127.1+0.5 is decidedly toward 
the outer Galaxy, so the X-ray sources are predominantly 
background galaxies. 

5. RESULTS 

Here we discuss the results of applying the techniques 
in § 4 to the first four SNRs from our Chandra sam- 
ple: SNRs G093.3+6.9, G315.4-2.3, G084.2+0.8, and 
G127.1+0.5. For each SNR, we first discuss its general 
properties, then details of the X-ray analysis, optical ob- 
servations and analysis, and finally identification of coun- 
terparts. 

5.1. SNR G093.3+6.9 

The radio source SNR G093.3+6.9 (also known as 
DA 530) was first identified as a SNR by Roger & Costain 
(1976), whose radio observations showed a shell 27' in di- 
ameter with bright rims and high polarization (Haslam, 
Salter, & Pauls 1980) — see Figure 2. The kinematic 
distance based on H I emission gives a distance of 2.5 ± 
0.4 kpc, but allowing for non-circular motion relaxes the 
distance limits and constrains the distance to be 1.0- 
3.5 kpc. 

While infrared and optical emission have been only 
marginally detected, SNR G093.3+6.9 has been detected 
by ROSAT (Landecker et al. 1999). Based on fits to the 
X-ray spectrum, Landecker et al. (1999) prefer a distance 
of 3.5 kpc. In X-rays, it appears superficially similar to 
SN 1006, which has been suggested to be the remnant of 
a Type la SN (Fesen et al. 1988, Allen, Petre, & Gotthelf 
2001). Taken together with the high Galactic latitude of 
7°, some authors believe that SNR G093.3+6.9 is also the 
young (« 5000 yr) remnant from a Type la SN (e.g., Lan- 
decker et al. 1999), but there are some problems with this, 
as none of the X-ray spectral models give fully consistent 
results and they imply that the SNR may have occurred in 
a especially low-density bubble such as might exist around 
a massive star (i.e. the progenitor of a Type II SN). 

Most of the X-ray fits in Landecker et al. (1999) assumed 
a hydrogen column density of TVjj = 2.1x 10 21 cm -2 , which 
is consistent with the value inferred from the H I data. 
This is the nominal value that we adopt for this source 
below. However, some of the fits required larger values of 
A^h, up to 6 x 10 21 cm~ 2 , while the total Galactic column 
density 14 in this direction is only 4 x 10 21 cm~ 2 . 

Lorimer, Lyne, & Camilo (1998) searched SNR G093.3+6.9 
unsuccessfully for a radio pulsar down to a limit of 
0.8 mjy at 606 MHz, implying a luminosity limit of 
Sd 2 < 9.8 mjy kpc 2 . Extrapolating to a frequency of 
1400 MHz using an average spectral index of —1.8 (Maron 
et al. 2000), this is still a factor of almost 10 brighter than 
the luminosities of the emerging class of faint radio pulsars 
such as PSR J0205+6449 in 3C 58 (Camilo et al. 2002b). 
Landecker et al. (1999) identified six point sources in the 
ROSAT PSPC data, three of which are inside the rem- 
nant (see § 5.1.1), but they conclude that none of these is 
likely the compact remnant of the explosion as they are all 
relatively hard (however, neutron stars such as the Crab 
can also have hard spectra). 

14 Determined using Colden, http:/ /asc. harvard.edu/toolkit/colden.jsp. 
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Table 17 

Classifications of X-ray Sources in SNRs G093.3+6.9, G315.4-2.3, G084.2+0.8, and G127.1+0.5 



SNR 


Stars 


Galaxies 


Uncertain a 


G093.3+6.9 


3 


15 





G315.4-2.3 


8 


4 


3: 17, 21, 28 


G084.2-0.8 


6 


8 


2: 5, 10 


G127.1+0.5 





9 


3: 12, 18, 22 


Total 


17 


36 


8 



Note. — Classifications follow the discussion in 
§§ 5.1.3, 5.2.3, 5.3.3, and 5.4.3. 

a Uncertain sources are those for which the clas- 
sification as a star or galaxy was unclear. This en- 
compasses both sources that definitely have coun- 
terparts but where the type is uncertain and 
sources where the counterpart itself is uncertain. 
The uncertain sources for each SNR are listed. 



5.1.1. X-ray Observations 

We observed SNR G093.3+6.9 with the Chandra X-ray 
Observatory (Weisskopf et al. 2000) on 2001 December 
17. Based on examination of the ROSAT/PSPC data, 
we determined the center of the remnant 15 to be at J2000 
position 20 h 52 m 14 s , +55°20'30", about 1.5' away from the 
nominal position from Green (2001). The aim-point was 
on the ACIS-I imaging array, and the final exposure time 
was 16.4 ksec. See Figure 2 for the placement of the ACIS- 
I detector relative to the SNR. A smoothed of the data is 
shown in Figure 6 with the X-ray sources that we identified 
(Tab. 4) labeled. 

The analysis of the Chandra data proceeded according 
to Section 4.2. In Figure 10 we plot the H-band counts 
vs. the L-band counts for the sources in Table 4. We 
also plot the lines for sources with power-laws having 
T = 0, 1.0, 1.5, 2.0 and blackbody models having kT^ — 
0.2, 0.4, 1.0 keV, all with A H = 2 x 10 21 cm" 2 , and the 
median and 25-/75-percentiles of count-ratios (Ch/Cl) for 
the sources from the CDF and Orion. Most of the sources 
are consistent with power-laws having indices from 0.0- 
2.0, such as one would expect for energetic pulsars like 
the Crab or for AGN. A few sources (SNR G093.3+6.9:l, 
SNR G093.3+6.9:7, and SNR G093.3+6.9:16) are softer, 
with implied temperatures (for blackbody models) of 
fcToo w 0.4 keV. 

Three of the sources were also detected by Landecker 
et al. (1999) in ROSAT data. Source SNR G093.3+6.9:l is 
source 6 from Landecker et al. (1999), SNR G093.3+6.9:2 
is source 5 from Landecker et al. (1999), and SNR G093.3+6.9: 
is source 4 from Landecker et al. (1999). While it is diffi- 
cult to compare directly because of the small numbers of 
counts involved, the count-rates and hardness ratios are 
roughly comparable between those that we observed here 

15 The geometric centers of SNRs were identified by eye, with the 
understanding that the actual site of the SNR is not always at the 
geometric center (e.g., Gvaramadzc 2002). 



and those predicted by converting the ROSAT count-rates 
to Chandra using W3PIMMS 16 

5.1.2. Optical/IR Observations 

We observed the field of SNR G093.3+6.9 a number of 
times with a number of instruments, as described in Ta- 
ble 9. The aim of the observations was to identify coun- 
terparts to the X-ray sources in Table 4 with progressively 
deeper exposures (§ 4.5). Data reduction for the optical 
data used standard IRAF routines to subtract the bias, flat- 
field, and then combine separate exposures. For the LFC 
and LRIS-B data, where significant focal-plane distortion 
prevented simple addition of data and where there are mul- 
tiple CCDs, we used the IRAF MSCRED package to flatten 
each image with custom distortion maps prior to addition. 
For the infrared data, we subtracted dark frames, then 
produced a sky frame for subtraction by taking a sliding 
box-car window of 4 exposures on either side of a reference 
exposure. We then added the exposures together, identi- 
fied all the stars, and produced masks for the stars that 
were used to improve the sky frames in a second round of 
sky subtraction. 

We determined BVRI photometric zero-points for the 
P60CCD data using observations of the Stetson fields 17 
PG1657 and NGC 6823 (Stetson 2000). We then deter- 
mined zero-points for the other optical data referenced to 
the P60CCD observations, employing appropriate trans- 
formations (Jorgensen 1994; Smith et al. 2002). 

We determined a K s zero-point for the P60IR data us- 
ing observations of the standard stars SJ 9101, SJ 9166, 
SJ 9177, and SJ 9188 (Persson et al. 1998). For the NIRC 
observations, we determined zero-points using 1-4 2MASS 
(Skrutskie et al. 1997) stars in each field (the K s mag- 
nitudes include a 0.3-mag systematic uncertainty arising 
from zero-point calibration). For the WIRC observations, 

16 See http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html. 

17 http://cadcwww.dao.nrc.ca/cadcbin/wdbi.cgi/astrocat/stetson/query 
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we determined zero-points using 400 unsaturated 2MASS 
stars. 

We performed absolute astrometry on the R-band 
P60CCD data, the LFC data, the LRIS data, and the 
WIRC data. After applying distortion corrections to the 
LFC and LRIS data (we did not use any distortion for 
the P60CCD data or WIRC), we solved for plate-scale, 
rotation, and central position relative to stars from ver- 
sion 2.2 of the Guide Star Catalog (GSC-2.2 18 ) for all but 
WIRC and relative to 2MASS stars for WIRC, getting 
residuals in each coordinate of 0722 (427 stars), 0708 (135 
stars), 0713 (613 stars), and O'.'IO (530 stars) for P60CCD, 
LRIS, LFC, and WIRC respectively. For the remaining 
data (ESI, P60IR, and NIRC) we used non-saturated stars 
from the other observations to boot-strap the astrometry, 
getting typical residuals of < 0705 in each coordinate. 

5.1.3. Notes on Individual Sources 

SNR G093.3+6.9:l This source is almost certainly an 
active star: the infrared colors are indicative of 
class GO or so, the source falls on the locus of stars 
in Figure 27, and the variability (Fig. 14) is typi- 
cal of that seen for stars (Marino, Micela, & Peres 
2000). The X-ray luminosity expected of such a 
star, 10 29-31 ergs s _1 (Katsova & Cherepashchuk 
2000) is implies that it is ~ 500 pc away, consistent 
with its K s magnitude (Fig. 32). 

SNR G093.3+6.9:2 This is likely a galaxy, given its 
position in Figure 27 and its constant X-ray flux 
(Fig. 14). 

SNR G093.3+6.9:4 This is likely a galaxy, given its po- 
sition in Figure 27. 

SNR G093.3+6.9:5 This is likely a galaxy, given its po- 
sition in Figure 27. 

SNR G093.3+6.9:6 The IR counterpart places it mid- 
way between the stars and galaxies in Figure 27. 
It could be a spurious counterpart, but there are 
other sources in its vicinity in Figure 27. Based on 
its hard spectrum and extremely red R — K s colors, 
it is likely a galaxy. 

SNR G093.3+6.9:7 While redder than SNR G093.3+6.9:l, 
it is otherwise similar, suggesting that this is very 
likely also an active star of type M2 or so. 

SNR G093.3+6.9:8 The optical counterpart is near the 
edge of the error circle, but it is far off-axis and 
therefore this is not unexpected. This source is oth- 
erwise consistent with being a galaxy. 

SNR G093.3+6.9:9 This source has the most extreme 
X-ray-to-IR flux ratio of the sources in Table 13, 
but it most likely is not the associated neutron star. 
The optical counterpart has a very small (072) offset 
from the X-ray source, suggesting that while it may 
be 2 mag fainter than other counterparts it is still 
likely to be a real association. It is still reasonably 
consistent with the CDF locus and it is somewhat 
red (R — K s > 3), contrary to known neutron stars 

18 http://www-gsss.stsci.edu/support/data_accoss.htm 



(Kaplan et al. 2003b) but similar to that of other 
galaxy candidates like SNR G093.3+6.9:10. Also, 
the X-ray spectrum is moderately hard, consistent 
with the power-law expected from an AGN. With- 
out an optical counterpart, it might be a candidate 
for a low-luminosity Crab-like pulsar (Lx would be 
w 2 x 10 31 ergs s _1 ), but with the counterpart it is 
very likely an AGN. 

SNR G093. 3+6.9:10 This is likely a galaxy, given its 
position in Figure 27. 

SNR G093. 3+6.9:13 This is likely a galaxy, given its 
position in Figure 27. 

SNR G093. 3+6.9:14 This is likely a galaxy, given its 
position in Figure 27 and hard spectrum. 

SNR G093. 3+6. 9:15 This is likely a galaxy, given its 
position in Figure 27. It may be extended, indicat- 
ing a low redshift (Fig. 17). 

SNR G093. 3+6.9:16 While redder than SNR G093.3+6. 9:1, 
it is otherwise similar, suggesting that this is very 
likely also an active star of type K5. 

SNR G093. 3+6.9:17 This is likely a galaxy, given its 
hard spectrum and position in Figure 27. The IR 
source may also be extended (Fig. 18). 

SNR G093. 3+6.9:19 There are multiple IR sources 
listed in Table 13, one just in the error circle and 
one outside (see Fig. 18). The true counterpart 
is likely the one in the circle (the first source in 
Tab. 13), which is likely a galaxy, given its position 
in Figure 27 and its colors. 

SNR G093. 3+6.9:24 This is likely a galaxy, given its 
position in Figure 27. 

SNR G093. 3+6.9:26 This is likely a galaxy, given its 
position in Figure 27. 

SNR G093. 3+6.9:30 This is likely a galaxy, given its 
hard spectrum and position in Figure 27, similar to 
SNR G093.3+6.9:19. The X-ray hardness ratio is, 
while harder than most galaxies, not unheard of in 
that context (Bauer et al. 2002). 

5.2. SNR G315. 4-2.3 

SNR G315.4-2.3 (RCW 86, or MSH 14-63) is a large 
(45'), non-thermal Galactic radio shell (Fig. 3) identified as 
a SNR by Hill (1967) and considered as the remnant of the 
historic supernova explosion SN 185 (Clark & Stephenson 
1977). 

The identification of SNR G315.4-2.3 with the historic 
supernova explosion of AD 185 is a matter of some con- 
troversy 19 . Geometric considerations would require a dis- 
tance to the remnant of less than 1 kpc (Strom 1994), 
while the kinematic distance of the remnant is found to be 
2.8 ± 0.4 kpc (Rosado et al. 1996). The latter distance 

19 While some recent reviews of the Chinese record even question 
the association of SN 185 with a supernova explosion event (Chin 
& Huang 1994; Schaefer 1995; Smith 1997), this is not universally 
accepted (Stephenson & Green 2002). 
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suggests a physical connection with an OB association 
(Westerlund 1969), an interpretation which is supported 
by the light elemental abundance which favors a Type II 
supernova and by the measured interstellar column den- 
sity. However, some recent modeling of X-ray data suggest 
a distance as close as 1.2 kpc and abundances more typ- 
ical of a Type la explosion (Bocchino et al. 2000), more 
in line with that expected if SNR G315.4— 2.3 were the 
remnant of SN 185, although these interpretations are by 
no means secure (Rho et al. 2002). In what follows, we 
assume the 2.8 kpc distance; if the SNR were closer (and 
younger), then our luminosity limits would be lower and 
even more constraining if SNR G315.4— 2.3 is the result of 
a core-collapse event. 

In addition to radio emission, SNR G315.4— 2.3 shows 
thin Balmer-dominatcd filaments (Long & Blair 1990; 
Smith 1997) and X-ray emission (Pisarski, Helfand, & 
Kahn 1984; Vink, Kaastra, & Bleeker 1997) that have the 
same general morphology. There is X-ray spectral varia- 
tion over the remnant, but the hydrogen column density 
is likely w 2 - 3 x 10 21 cm" 2 (Vink et al. 1997) (the total 
Galactic column density from Coldcn in this direction is 
9 x 10 21 cm" 2 ). The age of SNR G315.4-2.3 is somewhat 
uncertain, as it has features of both young (a few thou- 
sand years) and old (> 10 4 yrs) remnants (Dickel, Strom, 
& Milne 2001). On average, probably the best estimate 
for its age is 4000 yrs, although ages up to 10 4 yrs are 
not impossible (Rosado et al. 1996; Petruk 1999; Bocchino 
et al. 2000; Borkowski et al. 2001). 

Kaspi et al. (1996) searched SNR G315.4-2.3 unsuc- 
cessfully for a radio pulsar down to a limit of 1.3 mJy at 
436 MHz and 0.2 mJy at 1520 MHz, both implying a 1400- 
MHz luminosity limit of Sd 2 < 1.5 mJy kpc 2 . This is quite 
faint — a factor of three fainter than PSR J0205+6449. 

Based on ROSAT and Einstein data, Vink et al. (2000) 
identified an X-ray point source inside SNR G315.4— 2.3 
that they conclude is likely an active star. Unfortunately, 
this source is outside our X-ray observations so we can- 
not confirm or deny their conclusion. Similarly, Gvara- 
madzc & Vikhlinin (2003) used Chandra data of the bright 
south-west region to search for point sources, motivated by 
their hypothesis that SNR G315.4— 2.3 was a significantly 
off-center cavity SN. They find one source at the edge of 
the remnant (Fig. 3) without an optical counterpart, al- 
though the only limit they cite is from the Digital Sky 
Survey (mii m ~ 21), and conclude based on this and the 
source's X-ray spectrum that it may be a neutron star (we 
have examined 2MASS and see no clear counterpart to this 
source there, implying limits of J > 15.8 and K s > 14.6). 
However, these optical/IR limits are far from constraining 
(§ 5.2.2) and it is quite possible that, while not an active 
star like the other source they identify, this source is an 
active galaxy. In any case, all of these sources are quite 
far out from the nominal center of the remnant. If the ex- 
plosion occurred near what we define as the center of the 
SNR, none of the sources in Vink et al. (2000) or Gvara- 
madze & Vikhlinin (2003) could be the compact remnant 
without transverse velocities in excess of 1500 km s _ - 
not an unheard of velocity (§ 6.1.1), but certainly large. 

5.2.1. X-ray Observations 



We observed SNR G315.4-2.3 with Chandra on 2002- 
December-02. We determined the center of the remnant 
to be at (J2000) 14 h 42 m 50 s , -62°28'20". The aim-point 
was on the ACIS-I imaging array, and the final exposure 
time was 10.0 ksec. See Figure 3 for the placement of the 
ACIS-I detector relative to the SNR. A smoothed image 
of the data is shown in Figure 7 with the X-ray sources 
that we identified (Tab. 5) labeled. 

The analysis of the Chandra data proceeded ac- 
cording to Section 4.2. In Figure 11 we plot the H- 
band counts vs. the L-band counts for the sources 
in Table 5. Most of the sources are consistent with 
power-laws having indices from 0.0-2.0, but the sources 
are slightly softer than those in SNR G093.3+6.9. A 
few sources (SNR G315.4-2.3:2, SNR G315.4-2.3:6, 
SNR G315.4-2.3:8, SNR G315.4-2.3:9, SNR G315.4-2.3:26, 
and SNR G315. 4-2. 3:28) are soft, with implied tempera- 
tures (for blackbody models) of kT x < 0.4 keV. 

5.2.2. Optical/IR Observations 

For SNR G315.4-2.3 we also observed the field with a 
variety of instruments, listed in Table 10. The reduction 
proceeded as in § 5.1.2. 

We determined BVRI zero-points for the C40 data us- 
ing observations of the Stetson fields 20 E4 and L107 (Stet- 
son 2000). We then determined zero-points for the other 
optical data referenced to the C40 observations, employing 
appropriate transformations (Jorgensen 1994; Smith et al. 
2002). We determined a K s zero-point for the PANIC data 
with w 20 2MASS stars in each field. 

Astrometry was performed relative to 2MASS 21 . For the 
C40 and EMMI data, where there is some optical dis- 
tortion, we computed solutions (plate-scale, rotation, and 
central position) locally around each X-ray source, limit- 
ing the fields to ±1'. These solutions typically used 45 
stars and had residuals of 0'.'09 in each coordinate. For 
the MagIC and PANIC images we computed solutions for 
the entire image, using w 80 stars and giving residuals of 
0'.'05 in each coordinate. 

5.2.3. Notes on Individual Sources 

SNR G315.4-2.3:l This is likely a galaxy, given its po- 
sition in Figure 28 and its hard spectrum (Fig. 11). 
However, we only detect the counterpart in a single 
band and it is quite faint, so it is possible that the 
counterpart is a coincidence. 

SNR G315.4-2.3:2 This is likely a star, given its posi- 
tion in Figure 28 and its soft spectrum (Fig. 11). 
Based on its IR colors and magnitude, this may be 
a K2III giant star. 

SNR G315.4-2.3:3 The optical counterpart places it 
midway between the stars and galaxies in Figure 28. 
It could be a spurious counterpart, but there are 
other sources in its vicinity in Figure 28. Based on 
its hard spectrum, it is likely a galaxy. 

20 http://cadcwww.dao.nrc.ca/cadcbin/wdbi.cgi/astrocat/stetson/query 

21 The 2MASS data for this field were released before those for 
SNR G093.3+6.9, SNR G084.2-0.8, and SNR G127.1+0.5, where 
we used the GSC-2.2 for some of the astrometry. We have compared 
the results of 2MASS and GSC-2.2 astrometry and found them in- 
distinguishable. 
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SNR G315.4-2.3:5 This is likely a galaxy, given its po- 
sition in Figure 28 and its hard spectrum (Fig. 11). 

SNR G315.4-2.3:6 This is likely a star, given its po- 
sition in Figure 28 and its very soft spectrum 
(Fig. 11). Based on its IR colors, it is likely type 
KO or so. 

SNR G315.4-2.3:8 This is likely a star, given its posi- 
tion in Figure 28 and its soft spectrum (Fig. 11). 
Based on its IR colors, it is likely type G5 or so. 

SNR G315.4-2.3:9 This is likely a star, given its posi- 
tion in Figure 28 and its soft spectrum (Fig. 11). 
Based on its IR colors, it is likely type M7 or so. 

SNR G315.4-2.3:12 This is likely a galaxy, given 
its position in Figure 28 and its hard spectrum 
(Fig. 11). While the chance of a random star in the 
error circle is non negligible, only the reddest stars 
have R — K s > 6 (Fig. 32) and these stars make up 
only a small fraction of those in this line of sight 
(Nakajima et al. 2000). 

SNR G315.4-2.3:13 This is likely a star, given its po- 
sition in Figure 28. Based on its IR colors, it is 
likely type M6 or so. 

SNR G315.4-2.3:17 The optical/IR source is on the 
edge of the X-ray error circle in Figure 20 (to the 
west of center). There is also a faint smudge that is 
below the detection threshold and may be an exten- 
sion of the source just to the north of the error cir- 
cle. If the source to the west is the counterpart, the 
offset is large but not entirely unreasonable (Fig. 36; 
see also Barger et al. 2003). The source would then 
be a galaxy, given its position in Figure 28, hard 
spectrum (Fig. 11), and red colors (R — K s > 8; 
see Thompson et al. 1999). It is also possible that 
the northern source is the correct counterpart, in 
which case the source would also likely be a galaxy. 
A third possibility is that the X-ray and optical/IR 
emission are not entirely spatially coincident, such 
as would be the case for a nearby interacting binary 
star or a low-redshift galaxy. This might explain 
the X-ray-to-optical offset. It is unlikely that this 
source is a neutron star, as its spectrum is quite 
hard — harder than that expected of a Crab-like 
pulsar. 

SNR G315. 4-2. 3:19 This is likely a star, given its po- 
sition in Figure 28. Based on its IR colors, it is 
likely type F7 or so. 

SNR G315. 4-2. 3:21 There are two optical sources in 
the X-ray error circle (Fig. 21), or perhaps one ex- 
tended source. Separately, each source is near the 
star/galaxy boundary in Figure 28, and the X-ray 
spectrum is intermediate, so no quick classification 
is possible. This may be an interacting binary (CV, 
X-ray binary, etc.), sources that would not have 
been in the Orion sample but are present in the 
general Galactic/extragalactic populations (Hertz 
& Grindlay 1988; Grindlay et al. 2003). 



SNR G315.4-2.3:23 This is likely a star, given its po- 
sition in Figure 28. Based on its IR colors, it is 
likely type M5 or so. 

SNR G315.4-2.3:26 This is likely a star, given its po- 
sition in Figure 28. Based on its IR colors, it is 
likely type K5 or so. 

SNR G315. 4-2. 3:28 There are three possible counter- 
parts in Figure 21. The brightest is the closest to 
the center of the error circle and would be consis- 
tent with being a star, and it may be extended to 
the north, possibly indicating a binary companion. 
The source to the east of the circle's center is more 
consistent with a galaxy. Finally, the source to the 
south-west of the circle's center, on the edge of the 
circle, is also more consistent with a galaxy. Given 
the soft X-ray spectrum of this source, the bright 
stellar companion is likely the correct one. It is 
also possible that there is no companion (the ran- 
dom probability for a star is somewhat high, even 
for the closest source), in which case it would be a 
candidate neutron star, but a definitive answer will 
have to await additional X-ray observations (assess- 
ing the spectrum and variability and improving the 
position). 

5.3. SNR G084.2-0.8 

SNR G084.2-0.8 was identified as a Galactic SNR by 
Matthews et al. (1977), who observed a well-defined non- 
thermal radio shell (Fig. 4) with diameters of about 20' x 
14'. Feldt & Green (1993) identified CO and H I emission 
interacting with SNR G084.2— 0.8, giving the remnant a 
kinematic distance of 4.5 kpc and a size of 28 pc x 22 pc. It 
has not been detected in X-rays so there are no spectral fits 
to determine its age or temperature. The hydrogen column 
density to SNR G084.2-0.8 is iV H ~ 2 x 10 21 (again using 
Colden) and integrating the appropriate velocity range), 
while the total hydrogen through the Galaxy is Ah mix 
10 22 . Assuming a Sedov-phase remnant, we find (see § 4.1) 

1 /2 

U w 0.6T-. and since T7 w 1 holds for most SNRs, we 
can say that t A w 0.3-1.0 (with T 7 = 0.3-3). 

Lorimer et al. (1998) searched SNR G084.2-0.8 unsuc- 
cessfully for a radio pulsar down to a limit of 1.1 mJy 
at 606 MHz, implying a luminosity limit of Sd 2 < 
44 mJy kpc 2 . Extrapolating to a frequency of 1400 MHz 
using an average spectral index of —1.8, this is a fac- 
tor of almost > 50 brighter than the luminosity of 
PSR J0205+6449. 

5.3.1. X-ray Observations 

We observed SNR G084.2-0.8 with Chandra on 2002- 
November-24. We determined the center of the remnant 
to be at (J2000) 20 h 53 m 21 s , +43°26'55". The aim-point 
was on the ACIS-I imaging array, and the final exposure 
time was 24.6 ksec. See Figure 4 for the placement of the 
ACIS-I detector relative to the SNR. A smoothed image 
of the data is shown in Figure 8 with the X-ray sources 
that we identified (Tab. 6) labeled. 

The analysis of the Chandra data proceeded according 
to Section 4.2. In Figure 12 we plot the H-band counts 
vs. the L-band counts for the sources in Table 6. About 
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half of the sources are consistent with power-laws hav- 
ing indices from 0.0-2.0, but a number of the sources are 
significantly harder. These are likely distant AGN that 
have had their soft photons heavily absorbed by Galactic 
gas — a source with F — 0.5 and TVh = 2 x 10 22 cm~ 2 
would have HRl h = 0.8, similar to the hardest sources. 
A few sources (SNR G084.2-0.8:l, SNR G084.2-0.8:4, 
SNR G084.2-0.8:14, SNR G084.2-0.8:16) are soft, with 
implied temperatures (for blackbody models) of kT^ < 
0.4 keV. 

5.3.2. Optical/IR Observations 

We observed the field of SNR G084.2-0.8 a number of 
times with a number of instruments, as described in Ta- 
ble 11. The reduction proceeded as in § 5.1.2. 

We determined the photometric zero-points for the LFC 
data by bootstrapping from VRI Palomar 60-inch ob- 
servations of the Stetson fields 22 L98, NGC 7654, and 
PG 0231 (Stetson 2000) and employing appropriate trans- 
formations (Jorgensen 1994; Smith et al. 2002) to LFC 
observations of these fields. For the WIRC observations, 
we determined zero-points using 1700 unsaturated 2MASS 
(Skrutskie et al. 1997) stars. For the NIRC observations, 
we determined zero-points using about 60 stars from the 
WIRC image. 

We performed absolute astrometry on the LFC data and 
the WIRC data. After applying distortion corrections to 
the LFC data (we did not use any distortion correction 
for WIRC), we solved for plate-scale, rotation, and cen- 
tral position relative to stars from version 2.2 of the Guide 
Star Catalog (GSC-2.2 23 ) for all but WIRC and relative 
to 2MASS stars for WIRC, getting residuals in each coor- 
dinate of 0'.'13 (600 stars), and 0'/19 (2500 stars) for LFC 
and WIRC respectively. For the NIRC data we used non- 
saturated stars from the WIRC to boot-strap the astrom- 
etry, getting typical residuals of < 0'.'06 in each coordinate 
with about 60 stars. 

5.3.3. Notes on Individual Sources 

SNR G084.2— 0.8:1 This source is almost certainly an 
active star: the infrared colors are indicative of class 
F6 or so, the X-ray spectrum is quite soft, and the 
source falls on the locus of stars in Figure 29. 

SNR G084.2-0.8:2 This is likely a galaxy, given its po- 
sition in Figure 29 and its hard X-ray spectrum. 

SNR G084.2-0.8:3 This is likely a galaxy, given its po- 
sition in Figure 29 and its hard X-ray spectrum. 

SNR G084.2-0.8:4 This source is probably an active 
star: the infrared colors are indicative of class M0 or 
so, the X-ray spectrum is quite soft, and the source 
falls on the locus of stars in Figure 29. There is an 
offset between the X-ray and IR positions, but this 
is consistent with the X-ray positional uncertainty 
and the number of chance stars of this brightness is 
quite low. 

SNR G084.2-0.8:5 This is probably a galaxy, given 
its position in Figure 29 and its hard X-ray spec- 

22 http://cadcwww.dao.nrc.ca/cadcbin/wdbi.cgi/astrocat/stetson/query 

23 http://www-gsss.stsci.edu/support/data_access.htm 



trum. However, the multiple possible IR counter- 
parts make a clear identification impossible. 

SNR G084.2-0.8:6 This is likely a galaxy, given its 
very hard X-ray spectrum. It does not appear in 
Figure 29 as there are no L-band counts, but its K s 
magnitude is similar to that of other galaxies such 
as SNR G084.2-0.8:2. 

SNR G084.2-0.8:7 This is likely a galaxy, given its po- 
sition in Figure 29 and its hard X-ray spectrum. 

SNR G084. 2-0. 8:10 This is likely a galaxy, given its 
position in Figure 29 and its hard X-ray spectrum. 
However, the multiple possible IR counterparts may 
a clear identification impossible. 

SNR G084. 2—0. 8:11 This source may be an active star, 
given how bright its IR counterpart is (it is possi- 
ble, but not likely, that the counterpart is due to 
chance). However, the X-ray spectrum is quite hard 
and the IR colors are far redder than those typical 
for stars (an M7 star would have J — K s 1.3, 
which would require Ay ss 10 mag to get to the 
observed value of J — K s ). Therefore, this source is 
likely an active galaxy, where intrinsic J — K s > 2 
is not that unusual (e.g., Franx et al. 2003) and a 
larger foreground Ah is expected. 

SNR G084.2-0.8:12 This source is probably an active 
star: the infrared colors are redder than typical for a 
main-sequence star of its magnitude, but it could be 
a distant M5 giant. In addition, the X-ray spectrum 
is quite soft, and the source falls on the locus of stars 
in Figure 29. 

SNR G084. 2-0. 8:14 This source is almost certainly an 
active star: the infrared colors are indicative of class 
GO, the X-ray spectrum is quite soft, and the source 
falls on the locus of stars in Figure 29. 

SNR G084.2-0.8:15 This is likely a galaxy, given its 
position in Figure 29 and its hard X-ray spectrum. 

SNR G084. 2-0. 8:16 This source is probably an active 
star: the infrared colors are indicative of class Ml 
(or possible a K giant), the X-ray spectrum is quite 
soft, and the source falls on the locus of stars in 
Figure 29. 

SNR G084. 2-0. 8:18 This may be a galaxy, given its 
hard X-ray spectrum and the faint IR counterparts. 
However, the error circle is large, and the IR coun- 
terparts could be spurious. Even so, it would still 
probably be extragalactic in origin as the X-ray 
source is too hard to be Galactic. For a power- 
law source with V = 0.7 (very hard for a neutron 
star, but plausible for an AGN) we would need 
Nh > 5 x 10 22 cm~ 2 to give the observed hardness 
ratio. This value is significantly higher than the 
column density expected to SNR G084.2— 0.8, and 
is in fact even somewhat higher than the Galactic 
column density in this direction (although the AGN 
spectrum could have T < 0.7). 
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SNR G084.2-0.8:22 This source is probably an active 
star: the infrared colors are redder than typical for a 
main-sequence star of its magnitude, but it could be 
a distant M5 giant. In addition, the X-ray spectrum 
is quite soft, and the source falls on the locus of stars 
in Figure 29. 

SNR G084.2-0.8:23 This is probably a galaxy, given 
its hard X-ray spectrum and the faint IR counter- 
parts. 

5.4. SNR G 127. 1+0.5 

SNR G127. 1+0.5 (R5) is a 45'-diameter radio shell 
(Fig. 5), also known as R5, which was identified as an 
SNR by Caswell (1977). It is remarkable in that it has 
a bright (396 mJy at 1.4 GHz), flat-spectrum radio point 
source (G127.11+0.54) near the center (Caswell 1977, Jon- 
cas, Roger, & Dewdney 1989), but optical observations 
(Kirshncr & Chevalier 1978, Spinrad, Stauffer, & Harlan 
1979) and H I absorption measurements (Pauls et al. 1982; 
Goss & van Gorkom 1984) instead favor an identification 
with a massive elliptical galaxy at w 100 Mpc distance, 
similar to M87. 

The distance is likely 1.2-1.3 kpc, if it is in fact asso- 
ciated with the open cluster NGC 559 (Pauls 1977). It 
has not been detected in X-rays so there are no spectral 
fits to determine its age or temperature. The hydrogen 
column density to SNR G127.1+0.5 is N u w 2 x 10 21 
(again using Colden) and integrating the appropriate ve- 
locity range), while the total hydrogen through the Galaxy 
is TVjj ~ 1 x 10 22 . Again assuming a Sedov-phase remnant, 
wc find (see § 4.1) t 4 « 0.6T 7 " 1/2 , or t 4 « 0.2-0.8 (with 
T 7 = 0.3-3). 

Lorimer et al. (1998) searched SNR G127.1+0.5 unsuc- 
cessfully for a radio pulsar down to a limit of 0.8 mJy 
at 606 MHz, implying a luminosity limit of Sd 2 < 
1.4 mJy kpc 2 . Extrapolating to a frequency of 1400 MHz 
using an average spectral index of —1.5, this is a fac- 
tor of almost > 500 brighter than the luminosity of 
PSR J0205+6449. 

5.4.1. X-ray Observations 

We observed SNR G127.1+0.5 with Chandra on 2002- 
September-14. We determined the center of the remnant 
to be at (J2000) 01 h 28 m 32 s , +63°06'34". The aim-point 
was on the ACIS-I imaging array, and the final exposure 
time was 19.5 ksec. See Figure 5 for the placement of the 
ACIS-I detector relative to the SNR. A smoothed image 
of the data is shown in Figure 9 with the X-ray sources 
that we identified (Tab. 7) labeled. 

The analysis of the Chandra data proceeded according 
to Section 4.2. In Figure 13 we plot the H-band counts 
vs. the L-band counts for the sources in Table 7. Virtually 
all of the sources are consistent with power-laws having 
indices from 0.0-2.0. It may be that we have underesti- 
mated the 7V H to SNR G127. 1+0.5, which would put the 
distribution of hardness ratios closer to what is seen in 
the other SNRs, but this does not affect our analysis of 
the sources. 

5.4.2. Optical/IR Observations 



We observed the field of SNR G127.1+0.5 a number of 
times with a number of instruments, as described in Ta- 
ble 12. The reduction proceeded as in § 5.1.2. 

We determined the photometric zero-points for the LFC 
data from VRI Palomar 60-inch observations of the Stet- 
son fields 24 L98, NGC 7654, and PG 0231 (Stetson 2000) 
and employing appropriate transformations (Jorgensen 
1994; Smith et al. 2002) to LFC observations of the SNR. 
For the WIRC observations, we determined zero-points us- 
ing 900 unsaturated 2MASS (Skrutskie et al. 1997) stars. 
For the NIRC observations, we determined zero-points us- 
ing about 20-30 stars from the WIRC image. 

We performed absolute astrometry on the LFC data 
and the WIRC data. After applying distortion corrections 
to the LFC data, we solved for plate-scale, rotation, and 
central position relative to stars from version 2.2 of the 
Guide Star Catalog (GSC-2.2 25 ) for LFC and relative to 
2MASS stars for WIRC, getting residuals in each coordi- 
nate of Of! 12 (1200 stars), and Of! 22 (1100 stars) for LFC 
and WIRC respectively. For the NIRC data we used non- 
saturated stars from the WIRC to boot-strap the astrom- 
etry, getting typical residuals of < Of! 07 in each coordinate 
with about 20 stars. 

5.4.3. Notes on Individual Sources 

SNR G127.1+0.5:l This source is coincident with the 
non-thermal radio source G127.11+0.54 and with 
its optical counterpart (Kirshner & Chevalier 1978; 
Spinrad et al. 1979). The X-ray spectrum is quite 
hard, fitting with the identification as an active 
galaxy. Interestingly, this is the only source for 
SNR G127. 1+0.5 to be in the stellar locus in Fig- 
ure 30, when we know it to be extragalactic. The 
optical counterpart is somewhat brighter than those 
of most galaxies relative to the X-ray flux (similar 
to SNR G084.2-0.8:ll). This could be due to the 
extended stellar emission seen from this source, or 
due to the orientation of the optical jets. 

SNR G127.1+0.5:3 This is probably a galaxy, given its 
hard X-ray spectrum and the faint IR counterpart. 

SNR G127.1+0.5:4 This is probably a galaxy, given its 
hard X-ray spectrum and the faint IR counterpart. 

SNR G127.1+0.5:5 This is probably a galaxy, given its 
hard X-ray spectrum and the faint IR counterpart. 

SNR G127.1+0.5:8 This is probably a galaxy, given its 
hard X-ray spectrum and the faint IR counterpart. 

SNR G127.1+0.5:12 This source is among the softest 
of the sources in SNR G127. 1+0.5 and its position 
in Figure 30 is intermediate between the stars and 
galaxies, so it could be an active star. However, 
the spectrum is not actually all that soft, and it is 
more likely to be a galaxy, given its faint, red IR 
counterpart. 

SNR G127.1+0.5:13 This is probably a galaxy, given 
its hard X-ray spectrum and the faint IR counter- 
part. 

24 http://cadcwww.dao.nrc.ca/cadcbin/wdbi.cgi/astrocat/stetson/query 

25 http: / / www-gsss.stsci.edu/support / data_access.htm 
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SNR G127.1+0.5:14 This is probably a galaxy, given 
its hard X-ray spectrum and the faint IR counter- 
part. 

SNR G127.1+0.5:15 This is probably a galaxy given 
its hard X-ray spectrum and the faint IR counter- 
part. 

SNR G127.1+0.5:17 This is probably a galaxy given 
its hard X-ray spectrum and the faint IR counter- 
part. 

SNR G127.1+0.5:18 This source is among the softest 
of the sources in SNR G127. 1+0.5 and its posi- 
tion in Figure 30 is intermediate between the stars 
and galaxies, so it could be an active star. Like 
SNR G127.1+0.5:12, though, the spectrum is still 
reasonably hard. The source could be a galaxy, 
given its X-ray spectrum and the faint, red IR coun- 
terpart. If it were a star, it would have to be a 
distant M giant in order to produce the observed 
J — K s value. It could also be a chance coincidence, 
given how crowded the region is in Figure 26. 

SNR G127.1+0.5:22 This could be a galaxy, given its 
hard X-ray spectrum and the faint IR counterparts. 
The second counterpart in Table 7 is almost cer- 
tainly a star, as suggested by the last column of 
Table 7. The first source may also be spurious, but 
the chances of a real association are reasonable. 

6. LIMITS ON CENTRAL SOURCES 

In Sections 5.1.3, 5.2.3, 5.3.3, and 5.4.3, we showed that 
almost all of the X-ray sources from Tables 4-7 can be 
identified either with foreground or background sources. 
Therefore, there does not appear to be any detected neu- 
tron star in SNRs G093.3+6.9, G315.4-2.3, G084.2+0.8, 
or G127. 1+0.5. There are a small number of cases where 
either the association or the type of source (star versus 
galaxy; see Tab. 17) is uncertain, either due to an op- 
tical/IR detection in only one band and/or a detection 
at a somewhat large distance from the X-ray source, but 
there are certainly no sources that scream out "I am a 
neutron star." If we accept this, we can then draw two 
limits to the flux of any compact central source: a con- 
servative limit (Limit I), and a loose limit (Limit II). 
The conservative limit will be the flux of the bright- 
est source for which the optical/IR counterpart is at all 
in doubt. For SNR G093.3+6.9, this would be source 
SNR G093.3+6.9:8: this counterpart is faint, has unknown 
colors, is somewhat far from the X-ray source, and is some- 
what soft. For SNR G315.4-2.3, this would be source 
SNR G315.4— 2.3:1: again this is faint, has unknown col- 
ors, and is in a crowded region. For SNR G084.2— 0.8, the 
limiting source would be SNR G084.2-0.8:5: the multiple 
IR counterparts make a firm association impossible, and 
the spectrum is not so hard as to eliminate the possibility 
of a Crab-like pulsar. Finally, for SNR G127.1+0.5, the 
limiting source is SNR G127. 1+0.5:18, where the crowding 
and uncertain classification make a firm association impos- 
sible, and again the spectrum is not so hard as to eliminate 
the possibility of a Crab-like pulsar. Here we have played 
devil's advocate, and questioned all of the associations in 



Section 5 as much as possible. We in fact believe that the 
associations are reasonably secure, but we cannot be cer- 
tain. The looser limits in each case come from the faintest 
sources in Tables 4-7, and assume that all of the associa- 
tions in Tables 13-16 are correct. We present the limits, 
along with approximate fluxes and luminosities for three 
different source models, in Table 18. In Figure 37, we 
plot the blackbody limits along with the luminosities of 
the sources in Table 2. The limits for SNR G127. 1+0.5 
are significantly below those of the other SNRs as it had 
never been observed in the X-rays before, so we did not 
know what the level of the diffuse background would be 
and therefore selected an exposure time that would guar- 
antee sufficient counts from a source above even the most 
pessimistic background. 

Below we discuss the implications of not detecting any 
sources in the contexts of different models for what the 
sources could be (§ 2). We also include discussions of limi- 
tations imposed by our observing strategy, specifically the 
limited field of view and the gaps in the ACIS-I detector. 

6.1. Instrumental Limitations 
6.1.1. Field of View 

While we did not detect any bright X-ray sources 
that were obviously compact objects, it may be that 
this was because the sources had extremely high ve- 
locities that carried them beyond the ACIS-I field- 
of-view. This would imply v± > 1600 km s _1 for 
both SNR G093.3+6.9 and SNR G315.4-2.3, v ± > 
1700 km s" 1 for SNR G084.2-0.8, and v ± > 700 km s" 1 
for SNR G127.1+0.5 (assuming a centered explosion). 
v± = 1600 km s _1 is higher than the velocities of 99% 
of the radio pulsar population, while v± = 700 km s^ 1 is 
higher than 90% of the population (Arzoumanian et al. 
2002). For one source, this might be acceptable, but for 
two or more sources the chances become too low (~ 10~ 7 
for all four sources), requiring another explanation. 

Gvaramadze & Vikhlinin (2003) believe that SNR G315.4 
was the result of an off-center explosion in a cavity created 
by the SN's moving progenitor, and have similar hypothe- 
ses about other SNRs (e.g., Bock & Gvaramadze 2002). 

6.1.2. Chip Gaps 

The ACIS-I detector has gaps between the four CCDs 
where the sensitivity falls to zero. During normal observa- 
tions some compensation is made for this by the dithering 
of the spacecraft. For our observations, we ended up with 
approximately 5% of the area having an effective exposure 
that was ~ 50% of the nominal exposure. The wavdetect 
program used exposure maps to account for this effect 
when detecting sources, so that sources located in the chip 
gaps can have fewer detected counts but end up with the 
same significance as a source in the middle of the chips, 
so we should not be missing sources due to the chip gaps. 
But there is a 5% chance of having a source in the gap 
region that is not detected for which the flux/luminosity 
limits should be a factor of 2 higher than those stated in 
Table 18. 

6.2. AXPs 

The properties of AXPs in Section 2.3 allow us to state 
that there are no such sources in the central 8' of SNRs 
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Table 18 

Limits on Central Sources in SNRs G093.3+6.9, G315.4-2.3, G084.2+0.8, and G127. 1+0.5 



Model 


Limit I 




Limit II 






Counts a 


Fx ab 


L x ac 


Counts a 


fx ab 


L x ac 


SNR G093.3+6.9: 














BB (fc^ = 0.25 keV) 


25 


9.1 


2.6 


11 


4.0 


1.1 


PL (T = 1.5) 


11 


20.7 


3.7 


11 


9.1 


1.6 


PL (T = 3.5) 


)) 


13.1 


7.3 


11 


5.8 


3.2 


SNR G315.4-2.3: 














BB (fcToc = 0.25 keV) 


45 


24.7 


4.5 


11 


6.0 


1.1 


PL (T = 1.5) 


11 


55.8 


6.4 


11 


13.6 


1.6 


PL (T = 3.5) 


11 


35.3 


12.4 


11 


8.6 


3.0 


SNR G084.2-0.8: 














BB (fcToc = 0.25 keV) 


20 


4.3 


2.0 


12 


2.5 


1.2 


PL (r = 1.5) 


11 


9.9 


2.9 


11 


5.9 


1.7 


PL (r = 3.5) 




6.3 


5.7 




3.8 


3.4 


SNR G127.1+0.5: 














BB (fcToc = 0.25 keV) 


20 


5.4 


0.1 


12 


3.2 


0.1 


PL (r = 1.5) 




12.5 


0.3 


11 


7.4 


0.1 


PL (T = 3.5) 


11 


7.9 


0.2 


11 


4.8 


0.1 



a In the 0.3-8.0 keV band. 

b X-ray flux xl0~ 15 ergs s _1 cm~ 2 . 

c X-ray luminosity xlO 31 ergs s _1 , corrected for absorption according to 
§§ 5.1-5.4, and using the distances in Table 1. 



Note. — Fluxes and luminosities were computed using W3PIMMS. See § 6 
for definition of limit types. 
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Fig. 37. — X-ray luminosities (0.5-2 keV) as a function of age for neutron stars in SNRs from Table 2. Sources whose emission 
is primarily thermal are indicated with plus symbols, those whose emission is primarily non-thermal are indicated with stars, 
and those with only limits are indicated with triangles. The sources that have X-ray PWNe, which are typically > 10 times the 
X-ray luminosity of the neutron stars themselves, are circled. We also plot the limits to blackbody emission from sources in SNRs 
G093.3+6.9 (red hatched region), G315.4-2.3 (green hatched region), G084.2+0.8 (blue cross-hatched region), and G127.1+0.5 
(gold hatched region). A 30% uncertainty in the distance has been added to the range of luminosities given in Table 18 (i.e. 
we have taken the Type I limits with a 30% larger distance and the Type II limits with a 30% smaller distance, to give the 
widest probable range of luminosities) , and the likely range of ages is also shown. The cooling curves are the lp proton superfluid 
models from Yakovlev et al. (2003) (solid lines, with mass as labeled) and the normal (i.e., non-superfluid) M = 1.35 Mq model 
(dot-dashed line), assuming blackbody spectra and Roa = 10 km. These curves are meant to be illustrative of general cooling 
trends, and should not be interpreted as detailed predictions. The horizontal lines show the luminosity produced by blackbodies 
with Roo — 10 km and logToo (K) as indicated. Faster cooling than the curves is possible, due either to the presence of exotic 
particles in the NS core or to the full onset of direct Urea cooling for a heavier NS (e.g., Yakovlev et al. 2002b). 
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G093.3+6.9, G315.4-2.3, G084.2+0.8, or G127.1+0.5: 
the luminosity limits for T = 3.5 in Table 18 are at most 
10 32 ergs s -1 , or two orders of magnitude below those of 
most AXPs and a factor of 10 less than the "quiescent" 
states of the possibly variable AXPs. This discrepancy 
cannot be solved by a slight change in distance or absorp- 
tion, and is therefore quite firm. 

6.3. Cooling Radio-Quiet Neutron Stars 

Using the standard cooling curve (modified Urea only) 
in Figure 37, we would estimate luminosities of (2- 
5)x 10 32 ergs s -1 for any compact sources. Obviously, 
our limits are below those values. Our limits are also 
below the luminosities of most of the sources in Ta- 
ble 2: only PSR B2334+61, CXO J061705.3+222127, 
PSR B1853+01, and RX J0007.0+7302 are comparable 
(while PSR B1853+01 is dominated by non-thermal X- 
ray emission, any cooling radiation would have to be be- 
low this level). However, most of these objects are older 
(> 30 kyr) than the SNRs considered here (< 10 kyr), and 
three of these sources have substantial X-ray (and radio) 
PWNe that would make them detectable in the absence 
of point-source emission. 

With our limits in Table 18, any thermal emission from 
sources in SNRs G093.3+6.9, G315.4-2.3, G084.2+0.8, or 
G127.1+0.5 would have to have T m < 8 x 10 5 K (or kT^ < 
70 eV, for Roo rts 10 km): see Figure 37. Even for radii of 
3 km (§ 2.4), the limits on kT^ are w 100 eV. These are 
lower than expected from standard cooling curves (Page 
1998, Slanc, Helfand, & Murray 2002, Kaminker et al. 
2002; Yakovlev et al. 2002b) and would require some exotic 
physics (pion cooling, direct Urea cooling, etc.) that may 
be related to a more massive neutron star (Yakovlev et al. 
2003). 

6.4. Radio Pulsars 

If we compare the luminosity limits for the T = 2.0 
model for SNRs G093.3+6.9, G315.4-2.3, G084.2+0.8, 
and G127. 1+0.5 with the luminosities of known radio pul- 
sars (Tab. 3), we see that our limits are below most but 
not all of luminosities of pulsars found in SNRs. Trans- 
lating our limits to limits on E (roughly, Lx,o.i-2.4kcV ~ 
0.5L X: o.3-8kcV and usin g ^x,o.i-2.4kcV ~ 10~ 3 £; Becker 
& Trumper 1997) we find E < 3 x 10 34 ergs s" 1 (for 
the type I limit) or E < 8 x 10 33 ergs s _1 (for the 
type II limit) for SNR G093.3+6.9, SNR G315.4-2.3, and 
SNR G084.2-0.8, and E < 2 x 10 33 ergs s" 1 (for the type I 
limit) or E < 3 x 10 32 ergs s _1 (for the type II limit) for 
SNR G127. 1+0.5. Comparing with the values of E given 
in Section 2.1, these limits are below those of traditional 
pulsars, but are compatible with the values for the low- 
E HBPSRs. Our limits on Lx arc consistent with the 
one \ow-E HBPSR that has detected X-ray emission - 
PSR J1718-37184 (McLaughlin et al. 2003) — but its 
age and distance are poorly known and its luminosity of 
<~ 10 30 ergs s -1 is roughly what is expected from standard 
cooling curves for a source with r « 30 kyr in contrast to 
the sources discussed here (§ 6.3). 

Since we cannot constrain the existence of a low--E pul- 
sar in any of these four SNRs, we can ask what its period 
might be. To do so we must assume that initial period is 



much less than the current period so that the characteris- 
tic age r = P/2 P is similar to the actual age and that the 
braking index has the constant value n = 3. We know that 
this is not always the case (e.g., Murray et al. 2002b; Migli- 
azzo et al. 2002; Lyne 2004), but it is the best guess that 
one can make. Under this assumption one finds P <~ 1.4 x 
10 23 VT^(tE)-°- 5 s and B ~ 3.2 x 10 i2 ^(T 2 E)-°- 5 G, 
where I = 10 45 /45 g cm 2 is the moment of inertia. With 
t w 4000-6000 yrs, we find P > 3 s for SNR G093.3+6.9, 
SNR G315.4-2.3, and SNR G084.2-0.8, and P > 9 s 
for SNR G127. 1+0.5. These periods are larger than 
those of most but not all radio pulsars (Young, Manch- 
ester, & Johnston 1999, Camilo et al. 2000; McLaugh- 
lin et al. 2003), and may be high enough to take any 
source in SNR G093.3+6.9 or SNR G315.4-2.3 beyond 
the radio "death line." The implied dipole magnetic 
fields are also high, > 2 x 10 14 G for SNR G093.3+6.9, 
SNR G315.4-2.3, and SNR G084.2-0.8 and > 5 x 10 15 G 
for SNR G127.1+0.5, similar to PSR J1814-1744 and 
PSR J1847— 0130. This may be indicative of a growing 
population of such objects: young, non-energetic long- 
period pulsars. The lack of detected pulsed radio emis- 
sion in any of these four SNRs (Lorimer et al. 1998; Kaspi 
et al. 1996) may be intrinsic (i.e., there is no radio emis- 
sion), it may be an orientation effect, or it may just be 
that the SNRs have not been searched deeply enough over 
enough of an area, as there is now a growing number of 
radio pulsars with luminosities (defined here as -Fradiorf 2 ) 
below 1 mJy kpc 2 (e.g., Camilo et al. 2002b, a), far lower 
than typical for radio pulsars. 

6.4.1. Pulsar Wind Nebulae 

Any PWNe in the SNRs discussed here are > 3000 yrs 
old and may have already interacted with the reverse 
shocks (§ 2.2), so their sizes and brightnesses would be 
hard to predict. We therefore examine limits on PWNe 
for a range of sizes (as in § 4.3). We also scale to a fiducial 
size of 1 pc rj 1' — we did not detect any sources with 
those sizes in our images except for the known thermal 
emission from RCW 86. 

To quantify this, we take the limits on extended sources 
from Section 4.3 and Table 8. We convert the count lim- 
its to luminosity limits using a photon index of T — 1.5, 
getting the limits in Table 8. These limits are below the lu- 
minosities of virtually all young PWNe detected in X-rays 
(Possenti et al. 2002), but are consistent with some older 
sources such as the Vela PWN, CTB 80, and W44 (Pavlov 
et al. 2001a). However, these PWNe all have significant 
non-thermal radio emission, emission that is not present 
in SNR G093.3+6.9, SNR G315.4-2.3, SNR G084.2-0.8, 
or SNR G127. 1+0.5 since they are all shell-type SNRs. 

6.5. Binary Systems 

As many as 50% of massive stars originate in binary 
systems, which presumably give rise to X-ray binaries and 
eventually millisecond pulsars. One might expect that we 
could see a binary system where the more massive star 
has gone supernova but the less massive has not evolved. 
In this case, we would see X-ray emission that appears 
to be (but is not physically) associated with an optically- 
detected star that might be hard to distinguish from the 
active stars that make up the majority of Galactic X-ray 
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sources (if the star did evolve and the companion were 
close enough, it would donate matter to the compact ob- 
ject and would appear as an X-ray binary and would have 
different properties). 

However, after only one supernova the binary system 
would have a small space velocity, < 100 km s _1 (Pfahl, 
Rappaport, & Podsiadlowski 2002a, Pfahl et al. 2002b). 
It would therefore be restricted to a smaller region than 
the full search, which accommodates velocities up to 
1500 km s _1 (§ 6.1.1). A velocity of 100 km s _1 is an an- 
gular offset of w 30" (for an average distance of 3 kpc and 
age of 5 kyr) — much smaller than the ACIS-I field-of- 
view. 

As seen in Tables 4-7, we have found only one source 
within a radius of 30" of the center: SNR G127.1+0.5:l, 
which is known to be extragalactic. In fact, none of the 
stellar sources (as determined in §§ 5.1.3, 5.2.3, 5.3.3, and 
5.4.3) appear to be at distances beyond 1 kpc (comparing 
with Fig. 32). This is likely a selection effect of our X-ray 
flux limits: for stellar luminosities of 10 29 -10 31 ergs s _1 , 
our flux limits translate into distance limits of <~ 1 kpc. 

If none of the stellar sources could be companions to a 
neutron star, we must ask if any of the sources identified as 
galaxies could in fact be stars, and their identifications as 
galaxies could be coincidence. As seen from Figures 32-34, 
a main-sequence star at the distances of SNRs G093.3+6.9, 
G315.4-2.3, G084.2+0.8, or G127.1+0.5 would have K s w 
15-19, depending on stellar type (giant stars at these dis- 
tance would mostly be too bright, with K < 10). All of 
the "galaxies" in Sections 5.1.3, 5.2.3, 5.3.3, and 5.4.3 have 
K > 17 except SNR G084.2-0.8:ll (and this is too red 
to be at the distance/reddening of SNR G084.2-0.8), so if 
they were main-sequence stars they would be type K5 or 
later. However, even if one of these sources is a companion 
to a neutron star, the Type I limits still apply, and any 
neutron star would be under-luminous (§ 6.3). 

6.6. Accreting Black Holes 

If the SNc produced black holes (BHs) and not neu- 
tron stars, the black holes themselves would be invisible. 
X-ray emission would only be detected if there were ma- 
terial accreting into the BHs. Models for such emission 
are not very well understood. Thermal emission might be 
expected to come from the inner portion of the accretion 
disk itself, with an area of several times 7ri?| (Chakrabarty 
et al. 2001), where Rs is the Schwarzschild radius of the 
BH. For a 10 Mq BH (with R s = 15 km), the area would 
be > 1000 km 2 , or much larger than the limits on ther- 
mal emission (§ 6.3) for temperatures > 100 eV. It is 
also possible that the X-ray emission arises from Comp- 
ton scattering in an optically thin corona over a thin 
disk or via optically thin brcmsstrahlung emission from a 
hot advection-dominated accretion flow (ADAF, as in the 
model of Narayan, Barret, & McClintock 1997), but these 
models are not well enough developed to provide useful 
constraints. 

6.7. Type la Supernovae 

Type la SNe are believe to result from carbon detona- 
tion/deflagration of a white dwarf that has been pushed 
beyond its mass limit through accretion from a compan- 
ion star. The resulting explosion completely disrupts the 



star, synthesizing nearly a solar mass of 56 Ni which ulti- 
mately decays to Fe. No compact core is left behind. The 
SNRs from such events thus form a subsamplc in which 
we do not expect to find an associated neutron star. The 
mean rate for Type la SNe is ~ 20 — 25% of that for core- 
collapse supernovae (Cappellaro, Evans, & Turatto 1999), 
so we expect less than 20% of the observable SNRs to be 
the result of such explosions. 

The ejecta produced in Type la events differs consid- 
erably than that from core-collapse SNe. The former are 
rich in Fe and Si while SNRs from core-collapse events are 
richer in O and Ne. For young SNRs, the X-ray spectra can 
be used to identify those of Type la origin (e.g., Hughes 
et al. 1995, 2003; Lewis et al. 2003). Once the X-ray emis- 
sion is dominated by swept-up circumstellar or interstellar 
matter, such discrimination is much more difficult. SNRs 
originating from core collapse are often found near the 
molecular clouds in which the progenitors formed, while 
those from Type la events are from stars with sufficiently 
long lifetimes that they can have traveled far from their 
birthplaces. Thus, for SNRs located in the near vicinity 
of active star-forming regions one can reasonably assume 
that they originated from massive stars. However, the ab- 
sence of nearby star formation or dense molecular material 
is not necessarily a direct indicator of a Type la progen- 
itor. There are three Galactic SNRs that are commonly 
thought of as the products of la events — Tycho, and 
SN 1006, and perhaps Kepler (Baade 1945; Fesen et al. 
1988; Allen et al. 2001) — all reasonably young "historic" 
remnants (Stephenson & Green 2002). 

If any of the SNRs discussed here were known to be the 
result of a Type la explosion (as suggested by some authors 
for both SNR G093.3+6.9 and SNR G315.4-2.3), then we 
would not expect to see a compact remnant. Observing a 
known la remnant would then be a good test case for our 
methodology: as we have not found any candidate com- 
pact remnant, it demonstrates that our method is not in- 
clined to find false positives. However, while knowing that 
a SNR is from a Type la explosion precludes the existence 
of compact objects, the lack of compact objects does not 
require a SNR to be from a Type la. Therefore, we cannot 
conclude that SNR G093.3+6.9 or SNR G315.4-2.3 is the 
result of a Type la explosion. 

7. CONCLUSIONS 

There are 45 known SNRs that are reliably within 5 
kpc of the Sun. Most of these SNRs are expected to con- 
tain central neutron stars: < 20% are expected to result 
from Type la SNe and thus not contain a central com- 
pact source, while ~ 20% (dependent on the stellar initial 
mass function, the limiting mass for black holes, and bi- 
nary evolution; Heger et al. 2003) are expected to host 
a central black hole that may not be easily identified as 
such. Thanks to the persistent efforts by astronomers over 
the past four decades, central sources have been detected 
in the X-ray and/or radio bands in 18 of these SNRs, and 
three have been identified as probable Type la SNe. In 
some cases, only a centrally located PWN is detected, but 
in those cases (i.e. IC 443) it is reasonably assumed that 
the PWN is powered by a central compact source. 

We have begun a program of searching for compact 
sources in the remaining 23 SNRs. The program has been 
motivated by the discovery of a point-like X-ray source 
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at the very center of the youngest known Galactic SNR, 
Cas A. The program has two observational components: 
imaging with Chandra or XMM in the X-ray band fol- 
lowed by ground-based optical/IR followup. The latter 
is essential because of the high incidence of of interlopers 
(foreground stars and background galaxies). Such a com- 
prehensive program is possible thanks to the astrometric 
accuracy of the X-ray missions combined with high sensi- 
tivity. 

In this paper we report on Chandra ACIS-I observa- 
tions of four shell remnants (G093.3+6.9, G315.4-2.3, 
G084.2+0.8, G127. 1+0.5). We undertook optical and IR 
observations of every X-ray source detected with more 
than 10 counts. For all detections, we found, within astro- 
metric errors, a counterpart in one or more bands. These 
counterparts were consistent with either foreground (stars) 
or background (galaxies) sources. In particular, the X-ray 
flux and/or the X-ray to optical (IR) ratio of the detected 
X-ray sources were not as extreme as all known classes 
of neutron stars: accreting neutron stars, radio pulsars, 
AXPs and SGRs. 

In § 6 we discuss reasons why standard neutron stars 
were not found in these SNRs: they could have fallen in 
the gap between the chips, they could have escaped our 
field of view due to a very high velocities, they could be 
undetectable black holes, or they could not exist owing to 
the SNRs being the results of Type la explosions. All of 
these scenarios are unlikely for a single source, and even 
more so for all four, but are technically possible. If, on 
the other hand, these scenarios do not apply, then four 
remnants contain neutron stars that are fainter than our 
X-ray detection limit (typically, Lx < 10 31 erg s _1 in the 
0.5-10 keV band). 

We now consider this last (and most interesting) possi- 
bility. In the absence of other forms of energy generation 
(accretion, rotation power, magnetic field decay) the min- 
imum X-ray flux one expects is set by the cooling of the 
neutron star. From Section 6.3 and Figure 37 we immedi- 
ately see that the central neutron stars in these four rem- 
nants must be cooler than those present for example in the 
similarly-aged Puppis A, PKS 1205-51/52 and RCW 103. 

Our knowledge of the physics of cooling is by no means 
firm. There is considerable debate among theorists as to 
which of the multitude of physical processes can signifi- 
cantly affect the cooling output and as to which of the 
physical parameters (mass, rotation rate, magnetic field) 
controls these processes (Yakovlev et al. 2002a) . Nonethe- 
less, there is agreement that more massive neutron stars 
(with their larger mean densities) cool more rapidly than 
those with smaller mass; this expectation is illustrated 
in Figure 37. Thus our upper limits can be made con- 
sistent with the cooling possibility provided the central 
neutron stars in these four SNRs are more massive than 
known cooling neutron stars. Indeed, the known exam- 
ples of radio-quiet objects could well result from a strong 
selection effect, namely the earlier X-ray observations by 
Einstein and ROSAT detected the warmer cooling neu- 
tron stars (ignoring the neutron stars detected because 
of non-thermal emission). The existing data may already 
hint at a parameter affecting cooling, as exemplified by 
PSR J0205+6449 (Slane et al. 2002) and the Vela pulsar 
(Pavlov et al. 2001b), but also possibly by PSR B1853+01 



(Petre, Kuntz, & Shelton 2002) and RX J0007.0+7302 
(Slane et al. 2004). 

Of course, we also do not see rotation-powered pulsars 
such as the majority of the objects in Table 2. Recent 
observations are finding pulsars with lower radio luminosi- 
ties and values of E than ever before, and our limits would 
only be consistent with these newer sources. One might 
then ask why we see neither a standard cooling neutron 
star nor a standard active pulsar, assuming that there is 
no intrinsic correlation between these properties. It is pos- 
sible that there truly are no neutron stars in these SNRs, 
allowing one to speculate wildly about what actually is 
there. 

In subsequent papers we plan to report the X-ray ob- 
servations and ground-based follow up of the remaining 
19 SNRs. If no central sources are identified then the hy- 
pothesis proposed here, namely that there is a parameter 
that determines cooling of neutron star will be strength- 
ened. The simplest (and most reasonable) suggestion is 
that this second parameter is the mass of neutron stars, 
so that cooling observations could be used to "weigh" iso- 
lated neutron stars (as discussed by Kaminker, Haensel, 
& Yakovlev 2001). This hypothesis is at odds with the 
strong clustering of binary neutron star masses (Thorsett 
& Chakrabarty 1999). We note, though, that mass de- 
termination is only possible for neutron stars in compact 
binary systems where significant interaction with the com- 
panion may have taken place. Thus it is possible that 
neutron stars resulting from single stars are systematically 
more massive than those which evolve in compact binary 
systems (and have gained mass through accretion). 
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